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I.  INFORMATION  SYSTEMS 


A.  Project  6151.  STATISTICAL  DATA  PROCESSING  AND  PATTERN  RECOGNITION 

Principal  Investigator;  T.  M.  Cover 
Staff:  R.  King,  A.  El  Gamal,  W.  Rogers,  K.  Jarret , 
J.  van  Campenhout,  M.  Salehi,  N-H.  Duan 


1 . Objective 

Hie  purpose  of  this  project  is  to  investigate  problems  of  in- 
formation processing  and  pattern  recognition. 


2 . Current  Status  of  Work 

a . The  Arbitrary  Relationship  between  Probability  of  Error 
and  Measurement  Subset  (J.  van  Campenhout) 

Consider  a two-class  hypothesis-testing  problem  in  which 
n is  the  set  of  possible  features.  Let  S £ denote  a feature  sub- 
set and  P (S)  denote  the  optimal  Bayes  risk  (probability  of  error) 
e 

Incurred  idien  the  measurements  S are  used  . 

This  research  considers  the  possible  values  of  P (S)  as 

e 

a function  of  the  feature  subset  S in  a two-class  problem  with  equal 

prior  probabilities.  All  orderings  of  Pg(S)  and  S c Q satisfying 

the  natural  monotonicity  condition  S'  c:  S ^ P^(S')  > Pg(S)  can  occur. 

No  other  restrictions  exist  on  P (S),  0 < P (S)  < 1/2;  therefore,  the 

e e — 

known  result  can  be  extended  from  the  achievability  of  orderings  to  the 
achievabllity  of  numerically  specified  sequences. 

During  the  course  of  the  proof,  a family  of  multivariate 
normal  distributions  capable  of  inducing  any  allowable  numerically  spec- 
ified misclassif ication  probabilities  is  observed.  A paper  has  been 
submitted  for  publication. 


b.  Markov  Channel  Codes  (K.  Jarret) 

The  conjecture  that  first-order  Markov  encoding  is  suffi- 
cient to  achieve  any  rate  pair  within  the  capacity  region  of  a memoryless 
multiple-access  channel  with  feedback  has  been  disproved  by  counterexam- 
ple . 
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The  most  general  possible  first-order  Markov  encoding, 
with  four  free  parameters,  was  constructed  for  the  binary  erasure  chan-  • | 

nel  (x  ,x  binary,  y = x +x  ).  An  upper  bound  on  H(Y)  was  con- 
structed  and  maximized  over  the  code  parameters.  Ihis  maximum  H was 
less  than  R + R„  for  known  achievable  rate  pairs.  First-order  Mar- 

kov  encoding,  therefore,  is  not  sufficiently  rich  to  achieve  capacity.  I 

c.  Multiple-Access  Channels  with  Noisy  Feedback  (R.  King)  I 

Channel  capacity  can  be  increased  vihen  noiseless  feedback 
is  added  to  a multiple-access  channel — a channel  with  two  or  more  trans-  ^ j 

mitters  and  one  receiver.  In  a general  communications  network,  however,  j 

feedback  may  not  be  noiseless  but  only  statistically  related  to  the  chan-  I | 

nel  output.  Noisy  feedback  can  be  superior  to  noiseless  feedback  if  the  j 

corrupted  feedback  is  not  a degraded  form  of  the  output.  Total  coopera-  ..j 

! 

tion  is  often  possible  between  two  transmitters  sending  Independent  mes- 

\ I 

sages  to  one  receiver  even  if  the  observed  feedback  is  not  a deterministic  ^ 

function  of  true  output. 

t 

d . The  Feedback  Capacity  of  Degraded  Broadcast  Channels  i 

(A.  El  Gamal) 

The  capacity  region  remains  unchanged  if  feedback  is  added 
to  the  physically  degraded  broadcast  channel.  This  is  consistent  with 
Shannon's  observation  of  the  discrete  memoryless  channel  with  feedback. 

The  result  has  been  established  both  for  discrete  memoryless  and  contin- 
uous-amplitude gausslan  channels. 

e . The  Capacity  of  a Class  of  Broadcast  Channels 
(A.  El  Gamal) 

The  capacity  region  has  been  established  for  a class  of 
broadcast  channels  p(y,z*x)  for  which  I(X;Y)  > I(X:Z)  for  all  input 
distributions  p(x).  This  region  resembles  the  capacity  region  of  the 

broadcast  channel  with  degraded  message  sets  and  extends  the  work  on  ' ] 

degraded  broadcast  channels  . I 
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f.  The  Resolution  of  an  Apparent  Paradox  (J.  Van  Campenhout) 

The  peaking  phenomenon  of  the  Bayes  recognition  accuracy 
of  pattern  classifiers  with  unknown  underlying  statistics  is  being  in- 
vestigated . "nils  effect  is  caused  by  improper  comparisons  of  statisti- 
cally incomparable  models.  A formalization  of  comparability  has  been 
introduced,  and  some  of  the  results  obtained  in  the  literature  are  being 
analyzed  in  this  context.  A paper  is  in  preparation. 


B.  Project  6240.  OPTIMAL  CODES  WITH  VARIABLE  WEIGHT  SYMBOLS 

principal  Investigator:  J.  Gill 
Staff;  N.  Cot 


1 . Objective 

The  objectives  of  this  project  are  to  determine  procedures  for 
producing  optimal  uniquely  decipherable  codes,  using  symbols  of  unequal 
costs,  and  to  study  the  prop)erties  of  optimal  codes. 


2.  Current  Status  of  Work 


We  have  demonstrated  that  variable-length  coding  schemes  can 
be  applied  to  several  fields  other  than  communications . These  include 

(a)  game  theory  where  they  are  associated  with  the  determi- 
nation of  optimal  gambling  strategies 

(b)  data  allocation  in  storage  and  data  retrieval  [1]  where 
their  application  is  based  on  the  relationship  between 
double -chained  trees  and  variable-length  codes  ; a file 
can  be  organized  as  a double-chained  tree  in  which  keys 
correspond  to  variable-length  code  words 


In  (b)  , an  optimal  tree  (or  code)  corresponds  to  a minimal  average  search 
time.  If  0 < b^  < bg  < . . . <b^  are  the  costs  of  the  symbols  of  a code 
alphabet  and  if  A is  determined  by 
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then  the  converse  to  the  source-coding  theorem  provides  a lower  bound  on 

the  average  cost  c of  a code  with  code-word  probabilities 

[P,  . • • • .P„  }. 

1 n 


where  H(p, ,...,p  ) is  the  entropy  (in  bits)  of  the  probability  distri- 
1 n 

button  . 

1 n 

We  have  established  an  upper  bound  on  the  cost  of  optimal  codes  , 


H(p p ) 

1 n 

log^  >.  + bj^  + 6 


where  & < b^  is  a quantity  calculated  from  the  costs  b^,...,b^.  Using 
this  upper  bound,  we  can  efficiently  determine  nearly  optimal  codes  for 
unequal  probabilities. 

For  equiprobable  code  words,  a procedure  [2]  for  producing  op- 
timal codes  has  been  analyzed  and  shown  to  require  O(log^  n)  steps  to 
produce  optimal  codes  with  n code  words,  and  the  cost  of  these  code 
words  was  determined  in  Ref.  3.  This  procedure  is  being  extended  to 
nearly  equiprobable  code  words . 

Re ferences 

1.  L.  Stanfel,  "Optimal  Trees  for  a Class  of  Information  Retrieval  Prob- 
lems," Info  . Stor  . Retr  . , 9,  1973,  pp . 43-59. 

2.  N.  Cot,  "The  Complexity  of  the  Variable-Length  'Encoding  Problem," 
Sixth  Southeastern  Conference  on  Combinatorics,  Graph  Theory,  and 
Computing,  1975. 

3.  N.  Cot,  "cost  of  Optimal  Prefix  Codes,"  Johns  Hopkins  Conference  on 
System  Sciences,  Apr  1977. 


C.  Project  6302.  PROBLEMS  IN  CONTROL 

Principal  Investigator;  G.  F.  Franklin 
Staff;  G.  Aral,  S.  Shah 


'i 
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Objective 


1 . 

Tills  study  Is  concerned  with  design  problems  In  automatic  con- 
trol systems,  especially  when  model  uncertainty  Is  present. 

2 . Current  Status  of  Work 

Among  various  techniques  for  order  reduction  In  linear  constant 
systems,  partial  realization  has  the  following  disadvantages. 

• Based  on  experimental  sample  data  , the  results  are  strongly 
dependent  on  the  sampling  rate. 

• Because  Markov-parameter  partial  realizations  are  high-fre- 
quency approximations , they  are  Inappropriate  for  control 
purposes . 

To  overcome  the  first  difficulty,  the  following  approach  Is 
suggested.  Instead  of  taking  the  first  2n  Markov  parameters  (the  se- 
quence of  impulse  response),  the  sample  points  should  cover  several  main 
time  constants  of  the  plant.  One  can  take  parameters  spaced  k periods 
apart  so  that 

2nkT  = 3t 

where 

n = order  of  realization 

T = sample  period  of  data 

T = main  time  constant  of  the  plant 

This  is  a partial  realization  of  a sampling  rate  of  kT  which  is  often 
too  slow. 

Starting  with  the  transfer  function  H|^T(s),  one  of  the  same 

degree  [H  (s)l  can  be  constructed  at  a sampling  rate  of  T so  that  the 

th  ^ th 

(Ik)  Markov  parameter  coincides  with  the  1 Markov  parameter  of  Hj^^Cs). 

Unfortunately,  the  result  is  not  unique.  Among  the  multiplicity  of  solu- 
tions, however,  we  suggest  taking  the  one  with  the  slowest  poles.  An  al- 
gorithm to  calculate  H,j.(s)  from  been  developed  consisting 

of 
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• solving  for  the  zeros  of  a polynomial 

th 

• calculating  k roots  of  the  zeros 

• solving  a set  of  n-1  linear  equations 

A corresponding  method  for  the  multivariable  system,  based  on 
a series  expansion  of  matrix  expressions,  gives  rise  to  the  difficulty 
of  convergence.  Currently,  the  multivariable  algorithm  converges  only 
if  the  eigenvalues  of  the  sparsely  sampled  realization  satisfy 

'A  - 1 j < 1 


In  Fig.  1.1,  the  impulse  responses  of  the  plant 

G(s)  = j 

( z - 0.5) 

are  compared  to  the  first-order  partial  realization  based  on  Markov  pa- 
rameters using  a sparsing  of  4 and  the  second-order  partial  realization 
using  a sparsing  of  2.  The  more  preferred  low-frequency  modeling  is 
achieved  through  the  time  moments  of  the  plant,  which  are  the  coeffi- 
cients of  s in  the  power-series  expansion  of  H(s)  about  s = 0, 


m.  = 
1 


(-1)^  d^H(s) 


ds 


s=0 


(t)  t^  dt  = H(-F)"^“^G 


In  discrete  systems,  the  same  expressions  are  often  carried  over  to  the 
z domain  as 


(-1)^  dSl(z) 

Ml  = 

dz 


z=() 


which  does  not  determine  a physically  meaningful  quantity. 
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Fig.  1.1.  IMPUI^IE  RESPONSES  VS  FIRST-  AND  SECOND-ORDER  REALIZATIONS. 
Arrows  indicate  sample  Instants  where  the  Impulse  responses  of  the 
plant  coincide  with  that  of  the  model  (first  order,  second  order,  or 
both  are  indicated  by  the  numbers  at  the  arrows). 


A one-to-one  physical  relationship  of  discrete  power-series 
expansion  coefficients  to  those  of  the  continuous  system  can  be  estab- 
lished if  the  expansion  is  performed  about  z = 1 which  corresp)onds  to 
s = 0 in  the  s-plane.  The  moments  of  a discrete  plant,  therefore,  are 
defined  as 


= 


(-T)^  d^H(z) 


i: 


dz 


i * 

= - ■ /,  n(n  + 1)  . . . (u  + i - 1)  h 

1 • H 


(nT)  = H($  - 1) 


-1-1 


z=l 
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Discrete  partial  realizations  based  on  the  above  definitions  have  the 
properties  of  its  continuous  counterpart  in  that  they  assume  the  same 
steady-state  characteristics  to  polynomial  excitations  as 

2n-l 

u(t)  = ^ 

i=0 

D.  Project  6418.  OPTICAL  COMPUTATION  BASED  ON  RESIDUE  ARITHMETIC 

Principal  Investigator:  J.  W.  Goodman 
Staff;  A.  Huang 


1.  Objective 

Since  1 October  1976,  this  project  has  been  reactivated.  Its 
new  goal  is  to  explore  the  feasibility  of  realizing  an  optical  computer 
based  on  residue  arithmetic. 

2 . Current  Status  of  Work 

As  indicated  in  the  last  semi-annual  report,  funding  from  the 
Air  Force  Office  of  Scientific  Research  was  received  1 February  1977  to 
support  this  work.  During  the  month  of  January,  attention  was  focused 
on  the  decomposition  of  spatial  maps  into  cascades  of  simpler  maps  to 
perform  residue  calculations.  This  decomposition  allows  a trade-off 
between  the  number  of  separate  maps  that  must  be  selectable  at  a given 
point  in  the  processor  and  the  number  of  sequential  maps  required  to 
accomplish  a given  mathematical  operation. 

The  project  was  terminated  on  31  January,  and  work  was  contin- 
ued under  the  new  APOSR  grant . 


E.  Project  6502,  DATA -COMPRESS ION  TECHNIQUES  AND  APPLICATIONS 

Principal  Investigator:  R.  M.  Gray 
Staff:  
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1 . Ob Ject tve 


The  objective  of  this  work  is  to  develop  nearly  optimal  data- 
compression  algorithms  having  moderate  complexity. 

2 , Current  Status  of  Work 

With  the  joint  support  of  JSEP  and  AFOSR  and  during  the  past 
two  years,  approaches  to  trellis  and  tree-encoding  data -compress ion  sys- 
tems for  linear  gauss ian-source  models  have  been  studied  via  simulation. 
Such  a system  consists  of  a digital  nonlinear  filter  as  decoder  and  a 
matched-trellis  (or  tree)  search  algorithm  such  as  a Viterbi  algorithm 
as  encoder . 

Because  the  decoding  filter  determines  the  encoder  structure, 
the  design  of  this  system  is  based  on  the  selection  of  a decoder.  A 
natural  (and  not  original)  choice  is  to  take  the  decoder  of  a good  data- 
compression  system  (such  as  DPCM)  and  to  replace  the  encoder  by  its 
matched-trellis  search.  An  alternative  new  approach  is  the  "fake  pro- 
cess" wherein  a filter  is  designed  so  that,  when  driven  by  a memoryless 
binary  process,  the  output  "looks  like”  the  original  source.  Two  ap- 
proaches were  developed  for  generating  the  fake  process — one  based  on 
the  Central  Limit  Theorem  and  the  other  on  inverse-distribution  coders. 

Extensive  simulations  of  these  systems  were  conducted  for  gaus- 
slan  memoryless,  autoregressive,  and  moving-average  sources,  mean-squared- 
error  distortion  measures,  and  low-rate  codes.  The  trellis  codes  outper- 
formed the  traditional  DPCM  technique  by  1 to  2 dB  for  sources  with  memory 
and  the  traditional  PCM  by  approximately  0.4  dB  for  memoryless  sources. 


F.  Project  6601.  INFORMATION  THEORY  AND  INFORNUTION  PROCESSING 

Principal  Investigator:  M.  E.  Heilman 
Staff:  R.  Kahn,  D.  Andelman 


1 . Objective 

The  objective  of  this  project  is  to  investigate  the  relation- 
ships between  information  theory.  Information  processing,  and  complexity 
theory. 
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2 , Current  Status  of  Work 

Improved  coding  schemes  have  been  developed  for  the  wiretap  | 

channel  with  feedback.  Satellitization  schemes,  similar  to  those  used 
in  coding  for  broadcast  channels,  generate  achievable  (R,d)  curves  as 
opposed  to  single  (R,d)  points.  The  convexity  of  these  curves  has  been 
demonstrated  numerically,  and  a proof  is  sought. 

Generalizations  of  the  wiretap  channel  with  feedback  have  been 
developed.  In  one,  for  example,  there  are  two  competing  receivers.  The  ^ 

transmitter  has  a message,  but  neither  receiver  knows  to  whom  it  is  in- 
tended . The  transmitter  sends  some  data  , asks  for  feedback , then  sends 
more  data  . The  goal  of  each  receiver  is  to  ensure  that  it  properly  de- 
codes all  messages  intended  for  it  and,  subject  to  that  constraint,  to  | 

learn  as  much  as  possible  about  messages  intended  for  the  other  receiver. 

If  the  initial  transmission  does  not  indicate  which  is  the  intended  re- 
ceiver, there  is  hope  that  both  receivers  will  provide  correct  feedback 

and  that  the  transmitter  can  improve  performance  over  that  of  the  wire-  ' 

tap  channel  with  feedback  (which  is  equivalent  to  this  model  if  the  in- 
tended receiver  is  identified  at  the  start).  Somewhat  surprisingly,  it 
appears  that  it  may  not  be  possible  to  im.prove  performance.  Reasons  for 
this  are  sought,  and  it  is  hoped  that  these  results  will  aid  in  better 
understanding  security-capacity  trade-offs  in  communications. 

G.  Project  6701.  THE  RECURSIVE-FILTER  ADAPTIVE  LINE  ENHANCER'^'  | 

Principal  Investigator:  B.  Widrow 
Staff:  P.  Thompson 

1 . Objective 

The  adaptive  line  enhancer  (ALE),  consisting  of  a digital  fil- 
ter adapted  to  produce  a minimum  mean-square-error  prediction  of  a future 
input,  has  proven  to  be  a useful  tool  for  detecting  narrowband  signals 
obscured  by  wideband  noise  and  for  estimating  the  spectral  characteristics 

^This  work  is  supported  partially  by  the  United  States  Department  of  En- 
ergy. 


10 


I 


of  signals.  Earlier  research  has  provided  a good  understanding  of  the 
behavior  of  the  ALE  when  the  adaptive  filter  is  constrained  to  an  all- 
zero tapped -delay-line  (TDL)  configuration.  It  is  widely  recognized 
that  a recursive  filter,  containing  both  poles  and  zeros,  can  perform 
better  with  a smaller  number  of  parameters  than  an  all-zero  filter; 
however,  procedures  for  adjusting  the  parameters  of  such  a filter  are 
not  well  established.  The  objective  of  this  research,  therefore,  is 
to  develop  and  analyze  appropriate  adaptation  algorithms  for  a recur- 
sive-filter adaptive  line  enhancer.  A key  to  achieving  this  objective 
will  be  to  constrain  the  filter  to  a configuration  appropriate  for  the 
line-enhancer  application  in  which  high  "q"  resonances  are  desired.  By 
exploiting  this  implicit  prior  knowledge  of  the  optimal  filter  config- 
uration, it  will  be  ijossible  to  avoid  many  of  the  difficulties  normally 
encountered  with  the  adaptation  of  recursive  filters. 


2 . Current  Status  of  Work 

A possible  configuration  for  a recursive-filter  adaptive  line 
enhancer  (Fig.  1.2a)  is  obtained  from  an  all-zero  ALE  by  attenuating  the 
input  and  closing  a feedback  path  around  the  adaptive  TDL  filter;  this 
is  based  on  the  rationale  that  the  enhanced -signal  output  constitutes  a 
less  noisy  source  of  signal  with  which  to  drive  the  TDL  filter  than  does 
the  input.  The  poles  and  zeros  of  the  resulting  recursive  filter  are 
adjusted  simultaneously  by  adapting  the  weights  of  the  TDL  filter  in 
much  the  same  way  that  they  are  adapted  in  an  all-zero  ALE.  The  select- 
able parameter  a appearing  in  the  input  and  feedback  paths  controls 
the  degree  to  which  the  filter  transfer  function  is  influenced  by  its 
poles.  With  (X  = 0,  the  feedback  path  is  open  and  the  poles  play  no 
role;  with  a = 1,  the  input  is  cut  off  leaving  only  the  natural  modes 
which  are  dominated  by  pole  locations.  With  a < 1,  the  adaptive  fil- 
ter can  produce  the  desired  high  "q"  resonances  at  the  frequencies  of 
narrowband  inputs . 

The  system  in  Fig.  1.2a  is  best  understood  in  terms  of  its 
input-to -error  transfer  function,  written  as 


H(z)  = 1 - 


(1  - a)  A(z) 
1 - a A(z) 


1 - A(z) 

1 - a A(z) 


(1.1) 
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a.  Block  diagram 


w = SIGNAL  FREQUENCY 
c«^n=NYQUlST  FREQUENCY 

/ ''n^O.3 


■Re(z) 


b.  Loci  of  convergent  zero  locations  for 
second-order  case  where  0!  = 0.90 


Fig.  1.2.  A RECURSIVE-FILTER  ADAPTIVE  LINE  ENHANCER. 


where 


./X  “1  “2  -n  /,  ox 

A(z)  = a,z  + a„z  + ...  + a z (1.2) 

1a  n 

Ideally,  the  error  output  should  contain  all  of  the  noise  and  none  of 
the  signal  because  the  enhanced -signal  output  would  then  consist  of 
signal  without  noise.  At  convergence,  therefore,  the  input-to-error 
transfer  function  H(z)  should  approach  zero  gain  at  frequencies  where 
signals  occur  and  unit  gain  at  all  other  frequencies.  It  follows  that 
the  zeros  of  H(z)  should  approach  the  unit  circle  at  angles  correspond- 
ing to  signal  frequencies.  The  poles  of  H(z)  are  kept  "close"  to  the 
zeros  if  a = 1,  thereby  maintaining  a gain  of  approximately  1 at  fre- 
quencies away  from  those  of  the  signal. 

The  actual  convergent  zero  locations  computed  for  the  second- 
order  ALE  in  Fig.  1.2a  are  plotted  in  Fig.  1.2b.  The  curves  represent 
loci  of  zero  locations  for  varying  signal -to -noise  ratios  at  several  sig- 
nal frequencies.  It  can  be  seen  that,  for  sufficiently  high  SNRs , the 
zeros  approach  the  unit  circle  at  the  desired  frequencies;  at  lower  SNRs, 
they  converge  to  points  well  within  the  unit  circle.  It  is  also  apparent 
that  the  system  will  perform  much  better  at  signal  frequencies  near  half- 
Nyquist  than  at  frequencies  near  zero  or  (by  symmetry)  near-Nyquist  where 
the  zero  locations  tend  to  be  biased  away  from  the  signal  frequencies . 

For  larger  a,  the  zeros  would  move  closer  to  the  unit  circle  for  any 
given  SNR. 

An  alternative  configuration  for  a recursive  ALE  can  be  ob- 
tained by  extending  the  above  discussion  concerning  the  desired  form  of 
H(z).  It  can  be  argued  that,  in  addition  to  having  zeros  near  the  unit 
circle  at  angles  corresponding  to  signal  frequencies,  H(z)  should  have 
poles  at  the  same  angles  but  with  slightly  smaller  magnitudes  to  provide 
the  desired  sharp  nulls.  This  implies  that  the  input-to-error  transfer 
function  should  be  in  the  form  of 


H(z) 


1 - A(z) 

1 - A(z/p) 


1 


1 - 


-1  -2 

- a^z  - ^^7. 

— — 5 — 

pa^z  - p a^z 


- a z 
n 


n -n 

p a z 

n 


(1.3) 
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the  loci  of  convergent  zero  locations  are  plotted  for  a second-order  ALE 
with  this  configuration.  (We  use  p = 0.95  compared  to  a = 0.90  in 
Fig.  1.2b  to  provide  similar  pole/zero  separation  in  the  radial  direc- 
tion.) Figures  1.2b  and  1.3b  indicate  that  the  filter  configuration  as- 
sociated with  Eq.  (1.3)  will  place  its  zeros  more  accurately  at  low  sig- 
nal frequencies  and  low  SNRs  than  will  that  associated  with  Eq.  (1.1). 
Simulation  experiments  confirm  this  result. 

Both  recursive-filter  adaptive  line  enhancers  promise  to  pro- 
vide significant  signal  enhancement  and  accurate  spectral  estimates  with 
a small  number  of  adjustable  parameters.  Both  have  been  simulated  on  a 
minicomputer  and  were  found  to  converge  to  the  predicted  parameter  val- 
ues in  a reasonable  amount  of  time.  Future  work  will  be  aimed  at  quan- 
tifying behavior  and  comparing  performance  to  proven  signal-processing 
techniques.  Results  will  be  extended  to  higher  order  filters  and  multi- 
ple sinusoidal  signals  and  to  signals  with  nonzero  bandwidth  and  colored 
noise . 


H.  Project  6851.  IMAGING  OF  GAMMA-RAY  DISTRIBUTIONS 

Principal  Investigator:  A.  Macovski 
Staff:  A.  L.  Steinbach,  D.  Rosenfeld 


1 . Objective 

The  objectives  of  this  work  are  to  study  the  fundamental  char- 
acteristics of  a gamma-ray  source  distribution  and  to  derive  optimal 
methods  for  Imaging  such  distributions.  Because  refractive  effects  do 
not  exist  at  these  energies,  various  attenuative  masking  structures  are 
being  Investigated  with  respect  to  efficiency,  resolution,  depth  delin- 
eation, and  signal-to-noise  ratio. 

2 . Current  Status  of  Work 

Work  on  this  project  has  been  terminated.  Two  dissertations 
have  been  completed — each  describing  alternate  approaches  to  the  Imaging 
of  gamma-ray  distribution. 
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The  first  investigated  zone-plate  imaging  systems.  A complete 
deterministic  analysis  of  the  three-dimensional  Impulse  response  was 
made,  and  computer  simulations  served  as  the  basis  for  the  design  of  an 


optimal  configuration. 

The  second  focused  on  a modulated-aperture  imaging  system. 
Optimal  codes  were  derived  for  the  time  modulation  of  the  apertures,  and 
computer  simulations  were  made  to  illustrate  the  source-dependence  of 
these  codes.  The  stochastic  nature  of  the  imaging  process  was  analyzed 
in  detail . 
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I.  Project  7050.  STUDIES  IN  STATISTICAL  SYSTEM  THEORY 

Principal  Investigators:  T.  Kailath  and  M.  Morf 
Staff:  J.  Dobbins,  D.  Lee,  B.  Levy,  J.  Nickolls, 

G.  Verghese , E.  Verriest,  and  A.  Vieira 


1 . Objective 

The  purpose  of  this  project  is  to  study  the  basic  problems  in 
statistical  system  theory. 


2 . Current  Status  of  Work 

a . Realization  of  Time-Varying  Impulse  Responses 
(E.  Verriest,  T.  Kailath) 

Stationary  systems  have  many  properties,  and  the  theory 
has  been  developed  extensively.  The  known  results  obtained  from  station- 
ary systems,  therefore,  should  not  be  "lost"  in  the  time-varying  case, 
especially  if  the  system  is  "close"  to  stationary.  A measure  of  such 
closeness  is  defined  through  the  ( i/dt+3/ds) -opera tor , introduced  ear- 
lier in  our  study  of  the  distance  of  stochastic  processes  from  station- 
arlty.  Using  this  ojjerator,  the  time-varying  impulse  response  can  be 
expanded  in  a Taylor-like  series  around  stationary  functions.  Based  on 
a truncated  series,  a realization  of  an  impulse  response  is  obtained, 
consisting  of  time-invariant  systems  with  time-varying  coefficients  be- 
longing to  a certain  ring  of  functions.  Conditions  for  observability/ 
controllability  can  be  determined  in  terms  of  eigenvalues  of  the  invar- 
iant parts  . 
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Further  research  will  be  directed  toward  a more  abstract 
algebraic  theory  of  time-varying  systems — possibly  connecting  the  former 
results  with  those  obtained  by  Kamen . 

b . Unification  of  Scattering  Theory  and  Square-Root 
Approaches  to  Least-Squares  Estimation  (M.  Morf, 

T.  Kallath,  G.  Verghese , J.  Dobbins) 

Scattering  theory  has  proven  to  be  very  powerful  in  gen- 
erating such  new  least-squares  estimation  results  as  doubling  formulas, 
change-of-initlal-conditions  formulas,  and  backward  Markovian  models. 

The  square-root  approach  to  estimation  has  also  become  important.  A 
unifying  framework  to  relate  the  results  obtained  from  scattering  theory 
and  the  square-root  approach  would  be  valuable.  Not  only  could  the 
scattering  results  be  easily  translated  into  numerically  stable  square- 
root  algorithms  , but  the  combination  of  highly  visual  scattering-theory 
techniques  with  the  stochastic  Gram-Schmidt  processes  underlying  the 
square-root  approach  should  facilitate  the  development  of  insight  and 
new  results.  Some  progress  toward  unification  has  been  made. 

Change  of  initial  conditions  and  doubling  formulas  have 
been  derived  naturally  in  the  square-root  framework.  Translating  some 
scattering  results  to  square-root  form  apparently  requires  "square- 
roots"  of  indefinite  matrices.  Although  this  can  be  accomplished  through 
factorizations  in  the  form  of  M = A E A' , where  E is  a signature 
matrix,  stochastic  interpretation  of  these  factorizations  remains  in 
question  (it  could  lead  to  "imaginary"  random  variables). 


c . New  Algebraic  Methods  for  Two-Dimensional  Systems  Theory 

(B.  Levy,  S-Y.  Kung , M.  Morf,  T.  Kailath) 

To  extend  the  results  obtained  from  one-dimensional  sys- 
tems theory  and  one-dimensional  data  processing  to  higher  dimensional 
problems  (such  as  image  processing,  object  reconstruction,  and  geophys- 
ical and  meteorological  data  processing),  it  appears  that  new  mathemat- 
ical techniques  must  be  developed  . 

We  have  studied  the  properties  of  rational  matrix  transfer 
functions  depending  on  two  variables.  In  addition,  we  have  successfully 
extended  the  existing  one-dimensional  results  obtained  from  extracting 
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the  greatest  common  right  divisor  (GCRD)  of  two  2-D  polynomial  matrices, 
Sylvester  resultants,  and  matrix  fraction  description  (MFD)  of  2-D  ra- 
tional matrices.  We  have  used  these  results  to  define  state-space  con- 
trollability and  observability  in  a manner  consistent  with  the  minimality 
of  2-D  systems.  It  was  found  necessary  to  introduce  such  mathematical 
tools  as  "primitive  factorizations"  of  2-D  polynomial  matrices  and  some 
algebraic  geometric  concepts  that  appear  to  be  new  to  systems  theory  and 
specific  to  multidimensional  problems. 

We  have  also  developed  computational  techniques  for  the 
modeling  of  2-D  systems;  these  include  a generalization  of  Levinsons  al- 
gorithm for  the  recursive  estimation  of  a 2-D  stochastic  process  and  a 
generalization  of  Lanczo's  recursion  for  the  fast  Fade  approximation  of 
the  2-D  transfer  function,  given  its  Markov  parameters.  These  algorithms 
are  related  to  the  theory  of  2-D  or^^ogonal  polynomials  on  the  hyper- 
unit-circle or  on  the  hyper-real-line. 

Our  results  include  extensions  of  maximum-entropy  spec- 
tral estimation  techniques  to  the  2-D  case,  2-d  spectral  factorization 
methods,  and  new  smoothing  algorithms  for  2-D  stochastic  processes.  By 
noting  the  analogy  among  2-D,  integer,  and  differential-delay  systems 
(all  belonging  to  the  class  of  systems  over  a principal  ideal  domain 
ring),  we  were  able  to  obtain  new  realization  methods  and  matrix-frac- 
tion descriptions  for  such  systems. 

d . System  and  Restricted-System  Equivalence  of  Models 
Described  by  Polynomial  Operators  (T.  Kailath, 

M.  Morf,  S-Y.  Kung  , B.  Ije v y , G . Verghese) 

The  concept  of  system  equivalence  was  introduced  by  Ros- 
enbrock  to  generalize  the  notion  of  system  similarity  to  nonstate-space 
models  (described  by  polynomial  operators).  The  definition  of  equiva- 
lence is  relatively  simple;  given  two  descriptions  of  the  same  physical 
system  by  two  different  observers,  one  would  like  to  determine  whether 
these  descriptions  are  the  same  from  a "coordinate  free"  point  of  view. 

We  have  studied  equivalence  based  on  polynomial  operators 
anci  state-space  models.  Earlier  definitions  of  system  equivalence  have 
been  related,  and  an  alternative  definition  based  on  the  construction  of 
a complete  information  set  for  a given  differential  system  has  been 
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proposed.  The  canonical  decompiosition  of  the  state-space  model  Into 
controllable/observable  parts  can  also  be  extended  to  polynomial  sys- 
tems . 

In  some  cases,  however,  as  when  the  transfer  function  of 
the  system  is  not  correct  or  when  the  polynomial  matrices  are  not  in  a 
row  (or  column  reduced),  additional  restrictions  on  system  equivalence 
are  required  to  ensure  that  the  structure  oi  the  system  at  infinity  is 
preserved.  In  this  context,  it  appears  that  the  theory  of  2-D  homoge- 
neous polynomial  matrices  provides  a convenient  tool  for  the  study  of 
the  structure  of  jxilynomial  systems  at  infinity. 

e.  Model  Fitting  of  Covariance  Data  (A.  Vieira,  T.Kailath) 

If  a finite  segment  of  a covariance  function  is  known, 

one  can  fit  an  all-pole  (autoregressive)  model  to  it.  This  modeling  is 

useful  for  such  purposes  as  prediction  and  control.  It  can  be  realized 

by  solving  the  so-called  Yule-Walker  (or  normal)  equations.  Because  of 

their  Toeplitz  structure,  these  equations  can  be  solved  in  a recursive 

2 

and  efficient  manner  (n*"  opKjrations  for  an  n vn  matrix)  by  a Durbin- 
Levinson  algorithm  that  can  be  recognized  also  as  the  Szego  recursions 
for  polynomials  orthogonal  on  the  unit  circle.  We  have  established  the 
relationship  of  this  algorithm  to  some  synthesis  procedures  in  the  net- 
work theory  of  Brune  and  Darlington,  wherein  elementary  sections  of  a 
transfer  function  are  factored  out.  This  factorization  procedure  can 
be  applied  to  a class  of  transfer  functions  larger  than  all-pole  models. 
Based  on  this  knowledge,  we  were  able  to  develop  an  algorithm  for  the 
fitting  of  a pole-zero  model  to  covariance  data.  Details  are  still  be- 
ing Investigated  . 

f.  Ladder  Structure  for  Least-Squares  Estimation  (M.  Morf , 

A.  Vieira,  D.  Lee,  T.  Kailath) 

We  have  extended  our  earlier  ladder  structure  to  the  so- 
called  "nonwlndowed"  or  covariance  method  of  least-squares  estimation. 
Using  our  previous  results,  we  were  able  to  handle  the  case  where  data 
values  before  the  known  ones  are  considered  to  be  zero  (prewindowing). 
Although  this  has  produced  excellent  results,  windowing  of  the  data  is 
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often  considered  undesirable.  As  a result,  we  have  developed  a ladder 
structure  that  works  within  the  given  data  (no  windowing)  and  is  only 

slightly  more  complex  than  the  prewindowing  method.  A trade-off  between  ' 

complexity  and  better  results  must  be  decided  by  the  user. 

g . Extended  Interpretations  of  Ordered -Reduct ion  Algorithms 

(G  . Verghese , T.  Kailath) 

An  algorithm  developed  by  Dewilde  to  obtain  the  Smith- 
McMillan  structure  of  a rational  matrix  at  any  frequency  from  its  Laurent 
expansion  at  that  frequency  has  been  observed  to  be  identical  to  the 
"structure"  algorithm  of  Silverman.  The  structure  algorithm  is  used  in 
many  problems  in  system  theory  but  lacks  the  Smith-McMillan  interpreta- 
tion because,  in  most  of  these  problems  (such  as  system  inversion),  the  ^ 

frequency  of  interest  is  infinity  and  the  Smith-McMillan  form  applies 
only  to  finite  frequencies.  The  extension  to  infinite  frequencies,  how- 
ever, is  straightforward  and  necessary. 

/ 

Based  on  this  insight,  we  have  been  able  to  obtain  a uni- 
fied understanding  of  problems  as  diverse  as  immersion,  realization  from 
polar  expansions,  the  role  of  infinite-decoupling  zeros,  and  the  asymp- 
totic behavior  of  multivariable  root  loci  . , 

h.  A Scattering  Framework  for  Estimation  (B . Friedlander, 

T.  Kailath,  G.  Verghese) 

Our  earlier  study  of  a scattering  model  for  estimation 
focused  on  the  medium  properties.  Recently,  we  have  identified  the  v 

propagating  waves  in  the  medium,  which  are  the  smoothed  state  estimate 

and  an  ad  .joint  state.  As  a result,  we  now  have  simple  pictorial  ways  ■ 

of  visualizing,  deriving,  and  organizing  the  ma.jor  known  smoothing  for- 
mulas, Including  those  for  the  effect  of  altered  a-priori  data  on  the  ] 

estimates  . 

i.  Square-Root  Algorlthm.s  for  Parallel  Processing  (M.  Morf, 

J.  Dobbins,  T.  Kailath) 

1 

We  have  developed  a family  of  algorithms  to  process  data 
in  parallel  for  least-squares  estimation  problems.  In  a preliminary 
stage,  data  segments  from  dis.joint  time  intervals  can  be  processed 


20 


independently,  and  a compact  set  of  information  is  then  communicated  be- 
tween processors  to  yield  the  overall  optimal  estimates.  All  computations 
are  perfonned  in  square-root  form,  and  algorithms  have  been  develojjed  for 
both  filtering  and  smoothing.  In  the  preliminary  processing  stage,  fil- 
ters are  Independently  initialized  in  each  interval  so  that  Chandrasekhar- 
type  algorithms  offer  computational  reductions  even  if  the  models  are  only 
piecewise  constant.  Thus  far,  algorithms  have  been  developed  only  for  the 
known  model.  Work  is  still  in  progress  for  the  case  where  only  covariance 
information  concerning  the  observations  is  known. 


II.  DIGITAL  SYSTEMS 


A.  Project  6961.  DESCRIPTION  LANGUAGES  AND  DESIGN  FOR  GENERAL- 

PURPOSE  COMPUTER  ARCHITECTURES 

Principal  Investigators;  M.  J.  Flynn, 

W.  M.  vanCleemput 

Staff:  T.  Bennett,  D.  Hanson,  J.  Hupp,  R.  Lee, 

K.  Stevens,  S.  Wakefield,  J.  Iliffe, 

A.  von  Bechtolsheim 

1 . Objective 

The  purposes  of  this  study  are  to  develop  description  languages 
to  describe  computer  architectures  and  to  establish  a basis  for  under- 
standing the  limits  of  computer-architecture  design  (to  determine  the 
fastest  speed  for  the  execution  of  a program). 


2 . Current  Status  of  Work 

a . The  SPRINT  Printed -Circuit  Design  System 

(W.  M.  vanCleemput,  K.  Stevens,  T.  Bennett,  J.  Hupp) 

The  objective  of  this  system  is  to  develop  an  interactive 
computer-assisted  design  of  printed  circuit  boards.  The  SPRINT  system 
allows  for  the  manual  placement  of  critical  components  and  for  automatic 
placement  of  such  other  components  as  14-  and  16-pin  dual  in-line  pack- 
ages. The  interconnection  routing  module  manually  routes  critical  con- 
nections and  automatically  routes  noncritical  connections.  The  current 
system  is  limited  to  two  signal  layers;  however,  extension  to  multilayer 
boards  is  planned . 

The  input  to  SPRINT  is  called  the  Structural  Description 
Language  (SDL).  In  the  future,  the  SDL  will  also  be  used  as  input  to  a 
logic  simulator,  a fault  test-generation/simulation  system,  and  automatic 
logic  diagram  generation. 

SPRINT  is  implemented  in  MORTRAN  and  FORTRAN  IV  on  the 
IBM  370  at  SLAC  and  makes  use  of  the  Tektronix  4013  terminal.  The  SDL 
compiler  is  implemented  in  SPITBOL  and  SNOBOL  dialect.  The  output  is  in 
the  form  of  a plot  from  which  the  artwork  must  be  generated  manually. 


b . Computer-Aided  Layout  of  Large-Scale  Integrated  Circuits 

(W.  M.  vanCleemput , E.  Siutz,  W.  Marti) 

Although  several  systems  exist  for  the  automated  layout 
of  LSI  circuitry,  none  obtains  a layout  comparable  to  one  designed  man- 
ually. The  objective  of  this  project  is  to  develop  and  implement  a sys- 
tem that  uses  algorithmic  approaches  in  which  certain  decisions  are  the 
responsibility  of  the  designer.  It  is  expected  that  this  system  will 
reduce  design  time  considerably  at  no  expense  of  excessive  silicon  area. 
It  will  also  be  able  to  make  use  of  hierarchical  features  if  specified. 

c.  Interactive  System  for  Design  Capture  (W . M . vanCleemput , 

A.  von  Bechtolsheim , B.  Wong) 

In  the  current  design  system,  all  input  is  in  the  form  of 
the  SDL.  Because  designers  often  prefer  to  use  schematic  diagrams  to 
perform  their  tasks,  a system  is  being  designed  that  will  take  as  its 
input  a schematic  drawing  from  an  interactive  graphics  terminal  and  will 
output  an  SDL  description  for  further  processing  by  the  various  design- 
automation  programs . 

d.  Bounds  for  Maximal  Parallelism  (R.  Lee,  M.  J.  Flynn) 

This  is  a study  of  the  performance  limits  of  a single 
program  executed  on  a large  number  of  identical  processors  operating  in 
parallel  in  an  MIMD  (multiple  instruction  multiple  data)  organization. 
Eh^en  with  the  rapidly  decreasing  cost  of  LSI  microprocessors,  it  is  not 
economically  feasible  to  consider  vast  numbers  of  processors  within  the 
computer  architecture  to  speed  up  a computation  at  the  expense  of  com- 
ponent-processor efficiency.  Some  "acceptable"  level  of  efficiency  must 
be  obtained.  We  must  determine  the  overall  speedup  magnitude  that  can 
be  expected  by  increasing  the  number  of  processors  when  the  problems  of 
control  and  communication  that  accompany  the  cooperation  and  competition 
between  these  processors  are  ignored  . 

We  first  defined  a general  model  of  computation  on  a p- 
parallel  processor  by  distinguishing  between  logical  parallelism  (p» 
processors)  Inherent  in  a computation  and  physical  parallelism  (p  pro- 
cessors) available  in  the  computer  organization.  Wo  then  determined  such 
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performance  measures  as  execution  time,  speedup,  efficiency,  and  the 
space-time  product  so  as  to  evaluate  the  performance  improvements  (if 
any)  of  p-parallel  processor  systems  compared  to  uniprocessor  systems. 
The  results  obtained  indicated  that,  generally,  performance  depends  on 
both  computer  architecture  and  computation. 

We  then  derived  necessary  and  sufficient  conditions  for 
the  maximum  attainable  speedup  of  a p-parallel  processor  to  be  <mln 
(p/ln  p,  p*/ln  p*) . This  bound  has  never  before  been  analytically  es- 
tablished . 

In  addition,  if  there  are  always  sufficient  processors 
(P  > p*)  , conditions  c .n  be  derived  under  which  the  maximum  attainable 
speedup  of  a computation  is  p*/ln  p* , maximum  efficiency  is  1/ln  p* , 
minimum  execution  time  is  T^  X In  p*/p*,  and  the  minimum  space-time 
product  is  T^^  X In  p*  (where  T^^  is  the  execution  time  of  the  compu- 
tation of  a uniprocessor).  The  empirical  speedups  obtained  for  a large 
number  of  different  computations  indicated  that  80  percent  of  all  pro- 
grams examined  satisfied  these  conditions  and  had  maximum  speedups  of 
less  than  p*/ln  p* . 

e . Parallel  Information  Processing  in  Biological  Systems 

(S . Wakefield,  M.  J.  Flynn) 

The  interconnections  and  types  of  synapses  between  the 
units  of  a neural  subsystem  (such  as  the  stoma togastric  ganglion  of  the 
lobster)  have  been  determined  by  biologists,  as  have  the  stereotyped 
motor  patterns  that  it  produces;  however,  the  exact  mechanism  and  the 
sequence  and  duration  bounds  of  impulse  bursts  that  underlie  the  pro- 
duction of  the  coordinated  muscle-activation  patterns  are  unknown.  Such 
a mechanism  would  be  analogous  to  the  switching  network  responsible  for 
the  traversal  of  states  in  a digital  sequential  circuit.  Because  of 
this  analogy,  this  and  similar  biological  information-processing  sub- 
systems are  being  investigated.  In  addition,  speed,  power  requirements, 
size,  information  capacity,  and  other  characteristics  of  individual 
neurons  are  being  compared  to  those  of  electronic  information-processing 
components . 
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f.  Interpretive  Machines  (J.  K.  Illffe,  M.  J.  Flynn) 

The  motivation  of  this  work  Is  to  achieve  language  Imple- 
mentations that  are  more  effective  In  some  measure  of  translation,  exe- 
cution, or  response  to  the  user  than  could  otherwise  be  obtained.  The 
Implied  comjjarlson  Is  to  the  established  technique  of  compiling  a code 
Into  a fixed  general-purpose  machine  prior  to  execution.  Although  sub- 
stantial benefits  can  be  expected  from  microprogramming.  It  does  not 
represent  the  best  approach  to  design  when  such  contributing  factors  as 
wide  performance  range,  multiple  source  languages,  and  stringent  secur- 
ity requirements  must  be  satisfied.  An  alternative  Is  to  combine  Inter- 
pretation and  a primitive  Instruction  set  and  to  provide  security  at  the 
"microprogram"  level . 


B.  Project  7151.  COMPUTER  ARCHITECTURE 

Principal  Investigators:  E.  J.  McCluskey, 

S . S . Owlckl 

Staff:  D.  J.  Rossetti,  F.  W.  Terman 

1 . Objective 

The  objective  of  this  study  Is  to  gain  new  insight  Into  com- 
puter architecture  and  operating  systems.  We  intend  to  design  a small 
operating  system  and  to  verify  Its  correctness.  A second  goal  is  to 
develop  new  methods  to  evaluate  system  performance  and  to  use  the  re- 


sults obtained  for  better  system  design. 


2,  Current  Status  of  Work 


The  basic  methods  for  designing  and  verifying  the  correctness 
of  concurrent  programs  have  been  analyzed.  The  design  of  a probably 
correct  operating  system  Is  under  way,  and  the  high-level  description 
Is  nearing  completion. 

The  second  direction  of  this  project  is  focused  on  the  evalu- 
ation of  various  computer  architectures.  Two  methods  that  are  very 
helpful  In  measuring  performance  are  under  study — a highly  accurate 
modeling  technique  and  monitoring  systems . 
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A set  of  trace  programs  has  been  develop)ed  to  measure  micro- 
processor performance  and  to  produce  unique  information  concerning  the 
programming  ease  of  various  microprocessors . The  tracing  is  highly  pow- 
erful (it  can  record  interrupts  and  I/O  behavior)  and  almost  transparent 
(external  hardware  reduces  significantly  the  amount  of  perturbation  in- 
troduced by  the  tracer)  . 

The  investigation  of  interleaved  memory  is  being  extended  to 
the  effects  of  request  ordering  and  to  the  most  efficient  ordering,  A 
new  model  has  been  develojjed  to  evaluate  the  performance  of  gracefully 
degradable  multiprocessor  systems,  based  on  such  measures  as  availabil- 
ity, processing  power,  and  proportion  of  time  in  the  degraded  mode.  The 
best  trade-off  between  computing  power  and  availability  or  between  hard- 
ware vs  software  fault  detection  can  also  be  determined. 

a . Verifying  Concurrent  Programs  with  Shared-Data  Classes 

(S  . S.  Owicki) 

Monitors  are  a valuable  tool  for  organizing  operations  on 
shared  data  in  concurrent  programs.  The  mutually  exclusive  procedure 
calls  provided  by  monitors,  however,  are  sometimes  overly  restrictive. 
Such  applications  can  be  programmed  using  shared  classes  which  do  not 
enforce  mutual  exclusion.  In  one  method  for  verifying  parallel  programs 
containing  shared  classes  [ll,  one  proves  that  each  class  procedure  per- 
forms correctly  when  executed  by  itself  and  that  simultaneous  execution 
of  other  class  procedures  cannot  interfere  with  its  correct  operation. 
After  a class  has  been  verified,  calls  to  its  procedures  can  be  consid- 
ered uninterruptible  actions  ; this  simplifies  the  proof  of  higher  level 
program  comp)onents , Proof  rules  for  classes  and  procedure  calls  are 
presented  in  Hoare’s  axiomatic  style.  Several  examples  have  been  veri- 
fied, including  two  versions  of  the  readers  and  writers  problem  and  a 
dynamic  resource  allocator. 

b.  Computer  System  Measurement  (D.  J.  Rossetti) 

Performance  measurements  have  been  divided  into  three 
areas — large-scale  processor  instruction  tracing,  microprocessor  in- 
struction tracing,  and  characterization  of  large-scale  input/output 
environments . 
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The  large-scale  processor  instruction-tracing  facility 
is  complete  [2] . It  is  providing  data  for  various  studies  of  computer 


{jerformance , architecture,  and  operating  systems. 

Microprocessor  tracing  is  based  on  a hybrid  hardware/ 
software  scheme  to  collect  instruction  frequency-distribution  statistics 
for  microprocessor-application  programs.  Some  of  its  major  characteris- 
tics are  that  it 

• is  independent  of  the  program  being  traced 

• traces  bursts  of  instructions  at  a sample  rate  determined 
by  the  user 

• causes  minimal  perturbation  because  tracing  is  assisted  by 
the  hardware 

A third  area  of  work  is  the  gathering  of  sufficient  data 
to  characterize  the  input/output  environment  of  large  systems.  We  are  in 
the  process  of  obtaining  complete  I/O  trace  information  for  use  in  the 
architecture  of  novel  storage  systems  involving  new  addressing  schemes  , 
technologies,  and  system  structures.  Four  real-world  environments  are 
being  measured,  plus  a variety  of  operating  systems,  processors,  and  I/O 
configurations . 

c.  Memory  Interleaving  (F.  W.  Terman) 

A series  of  models  of  interleaved  memory  systems  has  been 
investigated  by  trace-driven  simulations.  The  basic  model  extends  the 
one  developed  by  Burnett  and  Coffman  [31  to  the  architecture  of  the  IBM 
360/370  with  its  variable-length  instructions  and  operands  [41.  It  also 
includes  multibyte  transfers  per  memory  access.  Variations  of  the  basic 
simulation  model  allow  the  study  of  the  effects  of  channel  interference, 
data-request  reordering,  and  data-queue  emptying  at  Interrupts.  The  sim- 
ulation model  has  been  extended  to  include  the  effects  of  more  than  one 
processor  accessing  the  same  memory. 

Memory  requests  for  the  simulations  are  obtained  from  two 
sets  of  instruction-by- instruction  trace  records.  The  first  set  traces 


370.  The  second  set  traces  samples  of  the  activity  of  the  CPU,  including 
supervisor  and  problem  states  . 

Preliminary  results  indicate  that,  for  high  degrees  of 
interleaving  and  large  memory  transfer  widths  (16-way  interleaving  and 
16  bytes  transferred  per  memory  access),  considerable  improvement  in  ef- 
fective memory  bandwidth  has  been  achieved  by  the  presence  of  the  addi- 
tional processors;  typical  values  are  a \'2  improvement  in  effective 
bandwidth  for  two  processors  and  a factor  of  2 improvement  for  four 
processors.  For  low  degrees  of  interleaving  and  small  memory  transfer 
widths  (four-way  interleaving  and  four  bytes  transferred  per  memory 
access),  the  additional  processors  reduce  the  effective  memory  bandwidth 


because  they  interfere  with  each  other  in  the  fetching  of  instructions. 
Fewer  instructions  are  fetched;  therefore,  fewer  instructions  are  decoded 
and  a smaller  number  of  operand  requests  are  generated  jjer  memory  cycle. 

The  simulation  studies  have  produced  some  interesting  in- 
formation concerning  IBM  360/370  programs.  For  example,  the  average 
instruction  length  is  3.5  bytes,  the  average  number  of  operand  bytes 
required  per  instruction  executed  is  2.5,  and  2/3  of  these  operand  bytes 
are  fetched  with  the  remaining  1/3  being  stored. 

Current  work  is  directed  toward  obtaining  a better  under- 
standing of  the  mechanisms  producing  the  effects  observed  from  the  vari- 
ous simulation  models.  A detailed  comparison  is  also  being  made  between 
the  results  from  the  simulation  models  and  the  various  theoretical  stud- 
ies in  the  literature. 
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III.  INTEGRATED  CIRCUITS 


A.  Project  5012.  IC  PROCESS  DESIGN  AND  COMPUTER  MULTI LAraR 

SEMICONDUCTOR  DEVICE  ANALYSIS 

Principal  Investigator:  R.  W.  Dutton 
Staff:  D.  B.  Estreich 

1 . Objectives 

Computer  methods  were  developed  to  predict  impurity  profiles 
after  processing,  and  physical  parameters  related  to  impurity  ciffusion 
were  calculated . 

Improved  models  will  be  developed  for  the  prediction  and  sim- 
ulation of  the  behavior  of  multilayer  semiconductor  devices  based  on 
layout  and  device-processing  history. 

2 . Current  Status  of  Work 

The  progress  of  integrated-circuit  design  has  been  toward 
higher  density.  The  merging  of  elementary  devices  to  form  new  semicon- 
ductor structures  has  been  one  method  to  achieve  higher  density.  A good 

2 

example  of  this  approach  is  integrated  injection  logic  (I  L)  vdiich  merges 
a lateral  PNP  with  an  inverse-operated  vertical  NPN  to  form  a functional 
logic  gate. 

Reducing  active  device  dimensions  is  another  method  for  in- 
creasing integrated-circuit  density.  In  some  cases,  undesired  parasltics 
exist  within  an  integrated-circuit  structure,  and  scaling  down  often  en- 
hances the  undesirable  features  of  these  devices.  The  best  example  may 
be  the  parasitic  NPN  and  lateral  PNP  transistors  inherent  in  standard 
CMOS  Integrated  circuits.  These  transistors  can  lead  to  a latch-up  fail- 
ure mode  (a  four-layer  p-n-p-n  structure  that  regenerates  into  a high- 
current  low-voltage  state).  A typical  solution  to  reducing  latch-up  has 
been  to  increase  the  lateral  dimensions  of  the  layout;  however,  this 
decreases  circuit  density. 

The  purpose  of  this  research  is  to  (1)  achieve  greater  under- 
standing of  the  physics  of  the  latch-up  failure  mode  and  to  quantify  the 
material  and  layout  parameters  that  control  latch-up  in  an  integrated 
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circuit  and  (2)  develop  techniques  and  procedures  to  analyze  and  predict 
IC  latch-up  susceptibility  from  knowledge  of  the  doping  levels  and  cir- 
cuit layout. 

An  important  aspect  of  this  research  is  to  develop  methods  to 
identify  potential  latch-up  paths.  This  requires  knowledge  of  the  cir- 
cuit layout,  including  bias  supply-contact  locations,  MOS  source  loca- 
tions, and  the  design  of  the  P-well  region.  Using  this  Information  plus 
substrate  resistivity  and  the  P-well  impurity  profile,  a computer  solu- 
tion to  Ftoisson's  equation  can  be  used  to  map  the  potential  distribution 
in  the  bulk  regions  so  as  to  identify  the  latch-up  path  (there  may  be 
more  than  one  in  some  layouts).  Measurement  procedures  are  also  being 
investigated  for  locating  the  latch-up  path,  using  thermal  methods. 

After  the  latch-up  path  has  been  identified,  computer  programs 
based  on  Gummel’s  algorithm  will  be  used  to  compute  the  current  gains  in 
the  parasitic  vertical  NPN  and  lateral  PNP  transistors.  This  requires 
information  concerning  the  impurity  profiles,  potential  across  the  base 
of  the  lateral  PNP,  and  bulk  minority-carrier  lifetimes.  From  the  cal- 
culated current  gains  of  the  two  parasitic  transistors,  latch-up  suscep- 
tibility can  be  predicted — the  critical  condition  for  the  initiation  of 
latch-up  is  that  the  product  of  the  NPN  and  PNP  current  gains  must  be- 
come unity. 

Spiecial  test  structures  have  been  developed  to  isolate  the 
geometrical  factors,  impurity  profiles,  and  carrier  lifetimes  and  to 
study  their  influence  on  latch-up. 

The  above  methods  and  computational  tools  should  be  applicable 
to  the  latch-up  analysis  of  any  Junction-isolated  IC.  This  work  should 
also  be  useful  in  considering  the  radiation  effects  in  CMOS  integrated 
circuits  . 
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IV.  SOLID  STATE 


A.  Project  5111.  TRANSPORT  PROPERTIES  OF  Al^Ga^.^As , AlxOa^.^Sb, 

AND  In  Ga,  As  SINGLE  CRYSTALS 
X 1-x 

Principal  Investigator:  G.  L.  Pearson 

Staff:  K.  Jew 

1 . Objective 

The  purpose  of  this  project  is  a systematic  study  of  the  trans- 
port of  free  carriers  in  these  ternary  systems . The  electronic  parameters 
are  free  carrier  mobility,  resistivity,  minority-carrier  d if fus ion  length , 
activation  energy,  and  trap  energy  levels.  The  distribution  coefficients 
of  selected  doping  impurities  as  a function  of  x are  significant  growth 
parameters.  The  results  obtained  should  prove  useful  in  the  design  of 
photolumlnescent  devices,  injection  lasers,  solar  cells,  infrared  detec- 
tors, and  bulk-effect  microwave  devices. 

2 . Current  Status  of  Work 

A series  of  (p)Al^Ga^  ^As  crystals  with  compositions  throughout 
the  0.0  < X < 1.0  range  was  grown  on  (lOO)-oriented  GaAs  substrates, 
based  on  liquid-phase  epitaxial  techniques.  Epitaxial  layers  grown  on 
semlinsulating  GaAs:Cr  substrates  were  used  for  Hall-effect  measurements. 
PN  junction  diodes  made  from  (p)Al^Ga j^_^As  layers  grown  on  (n^)GaAs;Te 
substrates  were  used  for  measurements  of  capacitance  voltage,  current 
voltage,  and  length  of  minority-carrier  diffusion. 

All  of  the  crystals  were  prepared  by  the  same  procedure.  Be- 
fore each  growth,  the  gallium  melt  was  baked  out  for  12  hr  in  flowing 
at  800 “C.  After  bakeout,  the  dopant,  GaAs  source,  A1 , and  GaAs  substrates 
were  loaded  into  the  growth  system.  The  p-type  dopant  was  Ge  having  an 
atomic  concentration  of  0.5  percent  in  the  melt.  The  amount  of  A1  was 
varied  to  produce  a range  of  ternary  compositions.  The  crystals  were 
grown  in  the  temperature  Interval  of  800°  to  780°C  and  at  a cooling  rate 
of  0.2°C/min.  The  crystal  compositions  of  the  Al^Ga^_^As  layers  were 
measured  through  electron-probe  microanalysis,  and  their  A1  distribution 
coefficient  and  growth  rate  were  obtained  . 
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Hall  measurements  of  the  (p)Al^Ga^  ^As:Ge  crystals  were  made 

using  the  Van  der  Pauw  technique.  Ciirrier  concentration,  mobility,  and 

resistivity  were  measured  at  77°  to  300 °K. 

Variation  of  the  hole  concentration  at  300°K  as  a function  of 

A1  composition  x is  plotted  in  Fig.  -1.1.  The  crystals  were  grown  from 

melts  containing  0.5  atomic  percent  Ge , and  the  effect  of  increasing  A1 

content  on  dopant  incorporation  can  be  observed.  The  carrier  concentra- 

18  -3 

tion  for  GaAs  at  x = 0.0  is  ^10  cm  \\tiich  decreases  exponentially 
with  increasing  A1  content  x as 

17 

p = 8.6  X 10  exp(-3.98x)  (4.1) 


ALUMINUM  FRACTION,  x 


Fig.  4.1.  CARRIER  VS  A1  CONCENTRATION 

X IN  LIQUID-PHASE  EPITAXIALLY  GROWN 

(p)Al  Gn,  As:Ge  SINGLE  CRYSTALS. 

.X  1-x 
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Because  of  the  rapid  decrease  in  hole  concentration  with  an  increase  in 
X,  samples  with  high  A1  content  are  closely  compensated  for  by  the  back- 
ground donor  concentration,  and  there  is  some  deviation  from  the  general 
behavior  described  in  Eq . (4.1). 

The  hole  mobility  at  300°K  as  a function  of  A1  composition  x 

2 2 

(Fig.  4.2)  decreases  from  170  cm  A^-sec  for  GaAs  to  40  cm  A^-sec  for  AlAs  . 
This  decrease  is  caused  by  the  increasing  hole  effective  mass  in  the 
Al^Ga^^  ^As  system  as  the  A1  fraction  x is  increased.  For  GaAs,  the  ef- 
fective mass  m*  = 0.48  m and,  for  AlAs,  m*=0.79  m . In  A1  Ga  As, 

p o pox  1-x 

it  is  assumed  that  the  effective  mass  varies  linearly  with  composition 
X so  that 

m*(x)  = (0.48  + 0.31x)  m (4.2) 

P o 
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Fig,  4.2.  HOLE  MOBILITY  VS  A1 
CONCENTRATION  x IN  (p)AlxGai_^ 
As:Ge  SINGLE  CRYSTALS  AT  300°C. 
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scatterinfj  is  dominant  at  room  temperature  and,  for  this  scat- 

-5/2 

terinp;  mechanism,  mobility  varies  as  m*  . Based  on  the  assumed  x 

P j 

dependence  of  m* , ' 

P 

U (x)  = A[(0.48  - 0.31x)  m (4.3) 

P o 


The  proportionality  constant  was  chosen  to  match  the  measured  mobility 

of  GaAs . Tlie  calculated  result  is  plotted  as  the  solid  curve  in  B’ig . 

1.2  and  is  in  good  agreement  with  experimentally  measured  values. 

Carrier  concentration  vs  temperature  generally  varies  as  p = 

p exp(-E  /kT) . where  the  magnitude  of  p decreases  with  increasing 
o a o 

.A1  composition  x,  as  described  above.  The  activation  energy  E (Fig. 
4.3)  becomes  greater  as  A1  composition  x is  increased  and  was  obtaineil 
by  curve  fitting  to  the  experimental  data.  In  the  0.0  <x  <0.6  range, 
activation  energy  increases  from  *=20  to  120  MeV  following  the  e.xpres- 
s ion  , 

E (x)  = 10  - 175x  (4.4) 

a 


Fig.  1.3.  ACTIVATION  ENERGY  OF  Ge 
IN  Al,,jGa^_,,j.As  V.S  ALCMI-NTM  COMPO- 
SITION X.' 
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Photoluminescence  spectra  of  these  crystals,  measured  in  the  direct  band- 


gap  range  (x  < 0.37)  yield  Ge  energy  levels  in  agreement  with  those 
obtained  from  Hall  measurements.  In  crystals  with  high  A1  content  (x  > 
0.6),  rhere  is  wide  deviation  from  Eq . (4.4)  as  a result  of  the  close 
compensation  between  the  Ge  acceptors  and  background  donors  so  that  car- 
rier concentration  vs  T does  not  vary  with  p = p^  exp(E^/kT) . 

Variation  of  mobility  with  temperature  is  similar  in  all  sam- 
ples throughout  the  A1  composition  range  (0.0  < x < 1.0);  however,  the 
magnitude  of  mobility  generally  decreases  as  the  A1  content  x increases. 

A maximum  u occurs  near  100 °K . Below  100 °K,  u decreases  because 
P P 

of  impurity  scattering;  above  100 °K,  it  decreases  as  a result  of  lattice 

-1 .3 

scattering  and  varies  with  temperature  as  T 

A series  of  (plAl^Ga^^  ^As:Ge  crystals  grown  on  (n^)GaAs:Te  sub- 
strates was  investigated.  Diodes  were  fabricated  from  As  grown  wafers, 
after  evaporating  the  metal  contacts,  by  scribing  and  cleaving.  Hie  ca- 
pacitance-voltage characteristics  of  these  hetero junctions  were  measured 
to  obtain  the  carrier  concentrations  in  the  crystal  and  in  the  built-in 
potential  between  the  two  hetero junction  materials  as  a function  of  com- 
position. The  hole  concentrations  obtained  by  the  C-V  technique  are  sim- 
ilar to  the  results  of  the  Hall  measurements  (Fig.  4.1).  The  current- 
voltage  relationships  of  the  diodes  were  used  to  analyze  the  transport 
of  carriers  across  the  (p)Al^Ga^  ^As-(n)GaAs  interface  as  composition  x 
was  varied . 

The  electron-diffusion  lengths  in  the  (p)Al^Gaj^  ^As  layers  were 

obtained  from  the  decay  of  the  short-circuit  current  which  was  induced  as 

the  electron  beam  of  an  SEM  was  scanned  across  the  cleaved  edge  of  the 

p-n  junction.  This  induced  current  decreases  exponentially  with  distance 

X as  I = I exp(-X/L  ),  and  decay  is  characterized  by  the  diffusion 
o n 

length  L . The  electron  diffusion  lengths  obtained  are  plotted  in  Fig. 
n 

4.4  as  a function  of  A1  composition  x.  At  x = 0,  = 6.3  pm  and  de- 

creases as  the  aluminum  content  is  increased.  Over  most  of  the  composi- 
tion range  x,  L « 1 uim. 

n 

Studies  of  (p)Al  Ga,  As  at  x > 0.5  are  being  continued.  This 
X 1-x 

work  will  then  be  extended  to  (n)Al  Ga . As:Sn  which  will  bo  characterized 

X 1-x 

via  Van  der  Pauw,  photoluminescence,  I-V,  C-V,  and  scanning  electi-on-micro- 
scope  techniques. 
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ALUMtNUM  FRACTION, X 

Fig.  4.4.  EIi:CTRON  DIFFUSION  LENGTH 

VS  A1  CONCENTRATION  x IN  (p)Alx 

Ga  As:Ge  SINGLE  CRYSTALS. 
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B.  Pro.iect  5244.  STUDIES  OF  HIGH-TRANSITION  SUPERCONDUCTORS  SUCH  AS 

V3Si  AND  N’b3Sn  ANT)  DEPTH  PROFILING  OF  Al^Gaj^_^As- 
GaAs  HETERO.IUNCTIONS  ^ 

I^lncipal  Investigator;  W.  E.  Soicer  , 

Staff:  I.  Linclau , J.  Miller,  C.  M.  Garner 

1 . Objective 

This  project  is  involved  in  two  areas  of  research.  The  first 
is  to  explore  the  electronic  structures  of  typo-II  (hard)  superconduct- 
ing alloys,  particularly  the  structure  of  Nb^Sn , to  determi.  e such  pa- 
rameters as  valence  bandwidth  which  is  important  to  the  theoretical  mod- 
els of  superconductors.  Emphasis  is  directed  toward  the  preparation  of 
atomically  clean  Nb  surfaces,  followed  by  evaporation  of  Sn  as  an  adlayer 
and  annealing  to  form  Nta^Sn . 

The  second  is  to  study  the  chemical  Interface  between  Al^Gaj^_^As 
and  GaAs  in  this  hetero.junction  system,  especially  the  chemical  width  of 

3H  ; 

1 

— ;J 


the  junction  and  surface  morphology  of  the  Interface  for  various  alumi- 
num concentrations  (0.3  < X < 0.85).  The  Schottky  barrier  Au-Af^Ga j^_^As 
interface  is  also  examined. 

2 . Current  Status  of  Work 

a.  Nb  and  the  Oxidation  of  Nb  (J.  Miller,  I.  Lindau) 

Because  some  of  the  highest  superconducting  materials 
are  Nb  and  Nb  compounds,  it  is  of  interest  to  study  the  relationship  of 
the  electronic  structure  of  these  materials  to  their  superconducting 
properties.  With  knowledge  of  the  valence-band  structure,  various  the- 
oretical predictions  associating  these  electronic  and  superconducting 
properties  can  also  be  tested.  In  addition,  the  oxidation  properties 
of  Nb  are  of  value  in  such  applications  as  microwave  cavities  and  Jos- 
ephson  junctions . 

In  our  experiment,  the  Nb  samples  were  heat  cleaned  in  an 
ultrahigh  vacuum  chamber  with  a base  pressure  of  1 X 10  torr . The  Nb 
was  heated  by  electron  bombardment  to  within  ~ 200 ° of  the  melting  point 
so  as  to  remove  all  contaminants.  Other  cleaning  techniques  have  been 
attempted  unsuccessfully.  Photoemission  was  applied  on  a He  resonance 
lamp,  using  hv  = 21.2  and  40.8  eV  lines  simultaneously.  Tlie  EDC  is  char- 
acterized by  two  d-band  peaks;  one  has  a large  density  of  states  near 

and  a shoulder  1.3  eV  below  E , and  the  other  appears  2.8  eV  below  E • 

F r 

The  measured  valence-band  width  (4  eV)  and  general  characteristics  of  the 

bands  are  in  agreement  with  the  bulk  density  of  states  calculated  by  Ho, 

Louie,  Chelikowsky,  and  Cohen. ^ 

Clean  Nb  was  then  exposed  to  small  amounts  of  oxygen,  and 

the  growth  of  oxide  peaks  was  monitored  via  UPS.  With  small  exposures, 

a peak  grows  at  6.5  eV  from  E which  is  related  to  a chemisorbed  stage. 

F 

At  larger  exposures,  a second  peak  appears  at  7.3  eV'  from  E and  con- 

I* 

tinues  to  Increase;  we  associate  this  with  the  formation  of  a surface 
NbO  (where  x = 2 is  a likely  candidate).  Higher  exposures  produce  a 

X 

third  peak  at  8.2  eV  which  grows  until  saturation;  at  this  point,  the 

^Phys.  Rev.,  B-15 , 1977,  p.  1755. 
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oxygen  peak  is  similar  to  a sample  left  in  air  and  known  to  form 

Saturation  exposure  is  associated  with  the  formation  of  Nb  O at  the 

^ o ^1 

surface.  Having  learned  to  prepare  clean  Nb , we  now  plan  to  study  the  ' 

Nb  alloys  because  of  their  high  superconducting  transition  temperatures. 

b . Auger  Profiling  of  III-V  Hetero, junctions  and  Schottky 
Barriers  (C.  M.  Garner,  C-Y  Su , W.  E.  Spicer) 

To  produce  III-V  hetero.junction  devices,  many  growth  j 

techniques  have  been  employed;  however,  the  reproducibility  of  these 
devices  was  found  to  be  poor.  Because,  in  most  III-V  hetero junctions  , 
the  energy  band  gap  will  change  in  going  from  one  side  of  the  interface 
to  the  other,  the  electronic  properties  of  a hetero junction  depend  on 

the  variation  of  composition  with  distance  through  the  interface.  Sim-  ^ 

ilarly,  Schottky-barrier  heights  often  vary  with  sample-preparation 
techniques.  To  study  these  problems,  we  have  used  ion-milling  Auger 
electron  spectroscopy  to  analyze  these  samples  for  variation  of  compo- 
sition with  depth. 

In  studies  of  Al^Ga^  ^As-GaAs  hetero junctions  grown  by 

liquid-phase  epitaxy  to  produce  the  sharpest  possible  interface,  it  was 

found  that  samples  grown  at  800®  and  750 °C  had  interface  widths  of  120 

and  100  A,  respectively.  On  the  other  hand,  samples  grown  by  LPE  had 

Interface  widths  of  approximately  500  to  1500  A;  however,  this  technique 

developed  fluctuations  of  A1  concentration  that  caused  a gradual  instead 

of  a sharp  increase  of  photoresponse  at  the  photon  energies  above  the 

band  gap  of  GaAs . Samples  of  vapor-phase  epitaxy  InP  grown  on  InGaAsP 

were  observed  to  have  large  concentrations  of  Ga  in  the  VPE  layer  of  InP, 

and  this  would  alter  the  performance  of  the  devices. 

In  studies  of  the  Schottky  barriers,  chemical  etching  of 

the  AlGaAs  samples  prior  to  deposition  of  Au  was  observed  to  increase 

the  thickness  of  the  o.xide  between  the  Au  and  A1  Ga,  As  and  to  change 

X 1-x 

the  barrier  height.  Anal/ses  of  Ag  on  heat-cleaned  InP  further  revealed 
that  a layer  of  In  appeared  on  the  surface  of  the  Ag  and  that  the  In  was 
in  the  form  of  an  In  oxide.  A much  smaller  amount  of  P was  also  found 
on  the  surface  of  the  Ag . This  would  drastically  affect  the  device  as 
a photoemitter  because  it  would  be  necessary  for  the  electrons  to  move 
through  this  In-oxide  layer  before  escaping. 
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Our  findings  indicate  that  growth  techniques  can  produce 
results  other  than  those  desired.  Additional  investigations  of  these 
systems  and  others  are  required  to  better  understand  the  growth  processes 
and  the  relationship  of  interface  structure  to  the  electronic  properties. 
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V.  RADIOSCIENCE 


A,  Project  3172.  DIGITAL  IMAGE  PROCESSING 

Principal  Investigator;  R.  N.  Bracewell 
Staff:  J.  Mandevllle 


1 . Objective 

The  purpose  of  this  project  was  to  Investigate  the  applications 
of  maximum-entropy  reconstruction  to  the  field  of  computerized  tomography. 
There  is  a consensus  that  the  benefits  of  iterative  techniques  are  con- 
fined to  circumstances  where  measurement  coverage  is  deficient,  for  exam- 
ple, ^^hen  projections  are  truncated  or  when  the  object  is  scanned  over  a 
limited  angular  extent.  In  situations  other  than  these,  a much  faster 
convolution  method  produces  images  of  acceptable  quality.  In  addition, 
the  trends  toward  faster  scanning  apparatus  and  larger  arrays  make  rela- 
tively slow  iterative  processing  undesirable  unless  the  processing  results 
in  substantial  improvements.  Because  of  these  faults,  effort  was  concen- 
trated on  comparison  of  an  error-tolerant  maximum-entropy  algorithm,  ART, 
and  the  convolution  method  in  deficient  measurement  situations. 

2 . Current  Status  of  Work 

Projection  data  sets  corresponding  to  full  coverage,  truncated 
projection  coverage,  and  limited  angular  coverage  have  been  computed,  and 
reconstructions  from  these  data  sets  were  then  computed.  The  four  recon- 
struction procedures  were 

• a maximum  entropy  algorithm 

• additive  ART 

• multiplicative  ART 

• convolution 

A major  conclusion  is  that  the  maximum-entropy  algorithm  does 
not  appreciably  advance  the  state  of  the  art  of  imaging  in  this  field  for 
the  following  reasons.  First,  it  has  been  proved  that  multiplicative  ART 
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converges  to  a maximum-entropy  image.  Not  only  does  this  relationship 
provide  a convergent  computer  algorithm  for  maximum-entropy  reconstruc- 
tion, but  it  also  suggests  that  the  many  empirical  tests  of  multiplica- 
tive ART  are  indicative  of  the  imaging  performance  of  the  maximum-entropy 
method.  The  reconstructions  computed  in  this  investigation  confirm  this 
qualitative  equivalence  between  multiplicative  ART  and  the  maximum-entropy 
a Igorlthm . 

Second,  it  appears  that  maximum-entropy  images  suffer  from  many 
of  the  same  artifacts  present  in  Images  computed  from  deficient  measure- 
ments by  simpler  and  faster  algorithms.  This  conclusion  is  confirmed  by 
empirical  tests  and  mathematical  analyses  of  the  equations  defining  max- 
imum-entropy images . 

Although  maximum-entropy  reconstruction  has  limited  utility  in 
computerized  tomography,  this  technique  offers  advantages  for  such  pro- 
cessing tasks  as  high-resolution  spectral  analysis  and  image  reconstruc- 
tion in  radio  astronomy.  We  have  continued  to  develop  the  iterative 
strategy  appropriate  to  these  uses.  Among  the  reasons  favoring  maximum- 
entropy  analysis  are  less  emphasis  on  processing  speed,  availability  of 
fast  algorithms  based  on  FFT , and  the  nonnegativity  constraint  has  more 
impact  on  data  reduction  than  it  does  in  computerized  tomography. 

B.  Project  3606.  RADIO  ACOUSTIC  SOUNDING  SYSTEM 

Principal  Investigator;  A.  M.  Peterson  i 

Staff;  M.  S Frankel , N.  Bhatnagar,  A.  Abuel-Haija  I 

1 . Objective  ; 

The  purpose  of  this  project  is  to  study  the  capabilities  and  1 

limitations  of  the  Radio  Acoustic  Sounding  System  (RASS)  as  a means  for 
measuring  temperature  profiles  and  winds  in  the  lower  atmosphere. 

2.  Current  Status  of  Work  i 

During  this  reporting  period,  important  theoretical  work  has 
been  completed,  which  delineates  the  effects  of  atmospheric  turbulence  j 

and  winds  on  the  performance  of  the  RASS  , Contrary  to  our  early  . 
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expectations  and  the  predictions  of  other  researchers  , our  study  has 

demonstrated  that  a radio-acoustic  system  can  operate  at  much  higher  | , 

frequencies  than  have  previously  been  employed  and  still  provide  data 
to  heights  of  greater  than  1 km.  The  significance  of  this  result  is 
that,  by  increasing  acoustic  frequency,  the  size  of  the  acoustic  source 
(as  well  as  the  RASS)  can  be  reduced.  With  a smaller  source,  it  will 
be  possible  to  build  a transportable  RASS  to  collect  real-time  atmo- 
spheric data  at  locations  of  the  user's  choice.  These  predictions  have 

been  experimentally  verified  [1]  at  an  acoustic  frequency  of  1 kHz  (440  i 

MHz  radio  frequency) . 

This  reporting  period  will  conclude  our  work  on  the  RASS  which 
was  initiated  as  a technique  for  real-time  temperature  profiling  the  | • 

first  few  kilometers  of  the  lower  troposphere  [2] . The  basic  physics 
of  the  RASS  is  as  follows. 

The  electromagnetic  refractive  index  of  air  (above  30  MHz)  in 
the  lower  troposphere  is  not  exactly  unity  but  is  a function  of  pressure  , ■, 

temperature,  and  humidity.  Because  of  its  pressure  dependence,  this  re- 
fractive index  can  be  altered  by  a short  pulse  of  sound  from  an  acoustic 
source.  An  electromagnetic  RF  signal  generated  by  a radar  passes  through 
the  sound  pulse  and  is  scattered  as  a result  of  the  Induced  refractive- 
index  variations.  This  scattered  signal  is  collected  by  a receiver  and 
is  processed  to  determine  atmospheric  temperature  as  a function  of  the 
acoustic-pulse  position  (range).  In  a static  atmosphere,  the  received 
signal  is  maximized  when  a Bragg  condition  is  established  between  the 
acoustic  and  electromagnetic  signals.  This  occurs  if  the  electromagnetic 
wavelength  is  twice  tl.e  acoustic  wavelength  and  results  in  an  in-phase 
addition  at  the  receiver  of  the  electromagnetic  signal  scattered  from 
successive  acoustic  wavefronts  in  the  pulse. 

A doppler  radar  measures  the  speed  of  the  sound  pulse.  Because 
the  speed  of  sound  in  air  depends  on  its  temperature,  the  temperature 
profile  of  the  lower  atmosphere  can  be  obtained  with  the  RASS  [31.  Fig- 
ure 5.1  illustrates  the  geometry  of  a bistatic  RASS  wherein  the  acoustic 
source  is  located  between  separated  RF  transmit  and  receive  antennas.  In 
a monostatic  RASS,  the  antennas  and  acoustic  source  are  located  at  the 
same  p)olnt . 
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Fig.  5.1.  REFLECTION  OF  ELEC- 


TROMAGNETIC ENERGY  FROM  AN 
ACOUSTIC  PULSE  FOR  A RASS 
B I STATIC  GEOMETRY. 

In  an  early  study  at  Stanford  University  [21,  the  reflection 
of  electromagnetic  energy  from  an  acoustic  pulse  propagating  vertically 
into  the  lower  atmosphere  was  analyzed  at  a carrier  frequency  of  85  Hz. 

A static  atmosphere  (no  turbulence  or  winds)  was  assumed,  and  no  attempt 
was  made  to  determine  the  feasibility  of  operating  the  RASS  at  acoustic 
frequencies  higher  than  85  Hz,  which  has  the  advantages  of  a smaller 
system  and  an  increase  in  the  resolution  of  the  measured  atmospheric 
parameters . 

Recent  work  considers  the  interaction  of  electromagnetic  and 
acoustic  waves  when  the  RASS  is  operated  in  a stochastic  environment 
characterized  by  turbulence,  winds,  and  mean-temperature  gradients. 

These  atmospheric  parameters  are  important  when  evaluating  RASS  perfor- 
mance because  they  could  affect  the  spatial  coherency  in  and  between 
the  acoustic  wavefronts  of  the  transmitted  acoustic  pulse.  This  decrease 
in  coherency  would  reduce  the  received  signal  levels. 

An  analysis  of  the  effects  of  turbulence  on  the  coherency  of 
the  acoustic  wave  reveals  that  the  acoustic-wave  amplitude  at  an  altitude 
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is  damped  by  exp[-(  |\l/j^(x,y^)  | )/2]  because  of  complex-phase  inter- 
ference. In  this  damping  factor,  first-order  pertur- 

bation in  the  complex  phase  of  the  acoustic  wave  caused  by  turbulence. 

These  perturbations  are  small  in  the  lower  troposphere  ; conse- 
quently, for  acoustic  pulses  with  carrier  frequencies  below  a few  kilo- 
hertz propagating  (under  typical  atmospheric  conditions)  to  altitudes 
of  a few  kilometers,  turbulence  has  little  effect  on  the  strength  of  the 
received  radio  signal.  If  altitude  x of  the  acoustic  pulse  is  much 

less  than  the  coherence  length  of  propagation  x , the  effect  of  turbu- 

c 

lence  on  RASS  performance  is  insignificant.  This  coherence  length  of 
propagation  is 
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\diere  m determines  the  acoustic  refractive-index  structure  parameter 

a 2 

in  the  presence  of  inhomogeneous  turbulence,  C is  the  refractive-ln- 

no 

dex  structure  parameter  at  altitude  x , k is  the  acoustic  wave  number, 

o a 

and  A is  a proportionality  constant  that  measures  the  outer  scale  of 

turbulence.  It  can  be  observed  that  the  coherence  length  varies  inversely 

with  frequency  and  the  strength  of  the  turbulence.  For  example,  consider 

a RASS  operating  at  an  acoustic  frequency  of  f^  = 2 kHz  and  assume  that 

= 10  ® m X =lm,  A=2,  m = 1.33,  and  the  velocity  of 

no  o Q y 

sound  c = 340  m/sec ; the  coherence  length  is  then  x = 4 x 10  m . 
o c 

This  value  implies  that  the  coherency  of  an  acoustic  wave  is  not  affected 
by  turbulence  throughout  the  \diole  atmosphere  when  f =2  kHz . 

Quiescent  atmospheric  conditions  (such  as  temjjerature  gradients) 
cause  the  dispersion  of  the  acoustic  wavetrain.  Under  these  circum- 
stances, a match  between  electromagnetic  and  acoustic  waves  can  be  ob- 
tained by  the  use  of  a modified  Bragg-scatter  condition.  This  match 
depends  on  the  temperature  difference  over  the  length  of  the  acoustic 
wavetrain:  however,  typical  temperature  gradients  in  the  lower  atmosphere 
(on  the  order  of  -6.5‘’K/km)  have  a negligible  effect  on  the  received 
power  of  a RASS.  The  effect  of  vertical  winds  on  the  strength  of  the 
received  radio  signal  can  also  be  negligible. 
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Mean  horizontal  winds  shift  the  focus  of  reflected  electromag- 
netic energy  from  its  origin,  resulting  in  a decrease  in  received  signal  | 

level  when  a monostatic  RF  system  is  used  . Because  the  acoustic  wave- 
fronts  act  as  large  spherical  reflectors,  the  principle  of  specular  re- 
flection, however,  can  be  utilized  in  a bistatic  radar  geometry  to  mea- 
sure this  wind  component  remotely  (in  addition  to  atmospheric  temperature  j 

profiles).  In  this  RASS  configuration,  the  radar  antennas  and  acoustic 

sources  are  aligned  in  the  wind  direction — the  transmit  antenna  upwind  j 

from  the  acoustic  source  and  the  receive  antenna  downwind.  This  technique 
has  the  potential  of  measuring  horizontal  wind  velocities  of  magnitude 
up  to  a few  tens  of  meters  per  second  in  the  lower  troposphere. 

An  experiment  fll  consistent  with  our  recent  theoretical  work  | 

was  conducted  at  an  acoustic  frequency  of  1 kHz  (440  MHz  RF) . For  peri- 
ods when  atmospheric  winds  were  less  than  =4  m/sec,  the  radar  antennas 
were  placed  diametrically  opposite  each  other  on  a circle  centered  about 
the  vertical  axes  of  the  acoustic  array.  A 6 m minimum  spacing  between 
these  antennas  w:is  chosen  so  that  the  RF  feedthrough  (leakage)  signal 
would  be  low  enough  to  avoid  saturation  of  the  RF  preamplifier.  Because 
this  separation  is  approximately  one  order  of  magnitude  less  than  the 
minimum  altitude  from  which  data  were  obtained,  this  configuration  is 
referred  to  as  a monostatic  RASS. 

When  winds  were  greater  than  =4  m/sec,  it  was  necessary  to 
configure  the  RASS  as  a bistatic  radar.  The  radar  antennas  are  aligned 
in  the  wind  direction — the  transmit  antenna  upwind  from  the  acoustic 
source,  and  the  receive  antenna  downwind.  This  configuration  of  anten- 
nas is  necessary  because  horizontal  winds  move  the  focus  of  the  acoustic 
wavefronts;  however,  by  placing  the  antennas  as  indicated  above,  a "spec- 
ular" reflection  can  be  obtained  when  the  focus  of  the  acoustic  wavefronts 
bisects  the  line  drawn  between  the  RF  antennas.  Under  specular  reflection 
conditions,  signals  comparable  in  magnitude  to  those  calculated  are  ob- 
tained when  the  acoustic  pulse  is  moved  in  the  windward  direction. 

Three  types  of  measurements  were  made  during  the  experiment. 

First,  data  were  collected  during  both  windy  and  quiet  conditions  to  de- 
termine the  effects  of  winds  or  the  maximum  range  of  the  system.  During 
quiet  conditions  (winds  less  than  a few  m/sec),  continuous  echoes  were 
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received  to  approximately  1 km,  using  a monostatic  RASS . When  winds 
exceeded  a few  m/sec,  radar  doppler  signals  from  the  acoustic  pulse  were 
limited  to  a 250  m altitude  by  wind  shear  near  the  ground.  By  orienting 
the  radar  antennas  in  the  direction  of  the  wind  and  moving  the  receive 
antenna  downwind,  however,  the  range  at  which  echoes  were  received  in- 
creased to  1 km.  In  addition  to  increasing  the  range  over  which  doppler 
signals  could  be  obtained,  the  bistatic  RASS  permitted  us  to  determine 
horizontal  wind  speed  and  direction.  These  data  were  obtained  by  measur- 
ing the  location  of  the  radar  antennas  and  the  delay  times  between  echoes 
received  at  each  antenna  location. 

Second , ground-level  temperature  measurements  were  made  with 
the  RASS  by  Bragg-matching  the  RF  and  acoustic  frequencies  at  low  alti- 
tudes. These  RASS  temperatures  were  compared  to  ground  measurements, 
and  good  agreement  was  observed . 

Third,  a measurement  was  made  to  understand  the  possible  degra- 
dation of  the  acoustic  wavefront  coherency  as  the  pulse  propagates  through 
regions  of  atmospheric  turbulence.  Results  support  our  theoretical  anal- 
ysis that  turbulence  does  not  significantly  degrade  the  performance  of  a 
high-frequency  RASS  over  at  least  the  first  300  m,  even  when  very  long 
acoustic  pulses  are  transmitted  . 
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1 . Objective 

The  objective  of  this  project  is  to  investigate  the  propagation 
of  ULF  electromagnetic  waves  in  the  sea  (frequencies  less  than  5 Hz)  from 
undersea  harmonic  dipole  sources,  taking  into  account  the  effects  caused 
by  the  presence  of  the  ocean  floor.  This  is  an  extension  of  earlier  work 
[1-31  dealing  only  with  an  ocean  of  effectively  Infinite  depth. 


2.  Current  Status  of  Work 


This  project  consists  of  (a)  derivation  of  the  mathematical  ex- 
pressions for  the  fields  of  harmonic  dipole  sources  located  in  an  ocean 
of  finite  depth,  (b)  evaluation  of  these  expressions,  and  (c)  analysis  of 
the  results.  The  progress  in  each  of  these  stages  is  reported  below. 


a.  Derivation  of  the  Expressions 

A complete  set  of  expressions  has  been  developed  to  deter- 
mine the  electromagnetic  fields  of  the  four  principal  elementary  dipole 
sources  (vertical  magnetic  and  electric  and  horizontal  magnetic  and  elec- 
tric) in  the  presence  of  a layered  conducting  medium.  Based  on  these 
formulas,  the  electromagnetic  fields  can  be  computed  in  any  layer  of  a 
multilayered  isotropic  medium  containing  any  of  the  above  dipole  sources. 
Both  displacement  and  conduction  currents  may  be  considered,  and  there 
are  no  restrictions  on  frequency  or  source-receiver  separation.  To  our 
knowledge,  a set  of  expressions  of  comparable  completeness  and  generality 
has  not  been  available. 

b . Evaluation 

These  expressions  contain  Sommerfeld  integrals  that,  gen- 
erally, cannot  be  evaluated  in  closed  form.  We  evaluated  the  integrals 
numerically,  therefore,  using  a technique  developed  in  our  Radioscience 
Laboratory  [41. 

c . Analysis  of  Results 

The  results  are  obtained  in  normalized  form,  with  distances 
expressed  in  units  of  skin  depth  & and  with  ocean-floor  conductivity 
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normalized  to  that  of  sea  water.  This  approach  is  economical  because 

the  frequency  parameter  is  absorbed  into  the  normalizing  factors,  and  j 

sea-water  conductivity  can  be  varied  if  desired  by  adjusting  these  fac- 
tors . 

To  obtain  a propagation  path  of  significant  length  between 
each  terminal  and  the  sea -a ir /sea-floor  interfaces  but  not  excessively 
attenuating  the  fields  in  the  water,  both  the  transmitter  and  receiver 
were  set  at  a depth  of  6/2  in  an  ocean  one  skin  depth  deep.  Represen- 
tative realistic  conductivity  contrast  ratios  (a„,  /a  of  10  ^ and 
2 floor  sea 

10  ) were  considered,  and  the  horizontal  source-receiver  separation  p 

-1  2 

ranged  from  10  to  10  6,  which  corresponds  to  separations  of  25  m to 

25  km  at  1 Hz . j 

From  our  analysis,  we  conclude  that,  for  p < 3 6,  the 
underwater  electromagnetic  fields  differ  less  than  3 dB  from  those  pro-  ; 

duced  by  a dipole  in  an  ocean  of  infinite  extent.  Thereafter,  the  com-  ] 

puted  fields  rapidly  become  very  much  greater  than  the  "infinite  ocean"  ^ ' 

fields  and  exceed  them  by  more  than  600  dB  at  p = 100  6.  This  indicates  | 

that  other  modes  of  signal  propagation  (such  as  the  up-over-and -down  mode  j 

above  the  surface  and  the  down-under-and-up  mode  below  the  water/floor  j 

I 

interface)  must  be  the  principal  mechanisms  that  carry  the  signal  to  these 
large  distances.  We  have  developed  a method  for  separating  the  contribu- 
tions from  these  modes  and  the  contribution  from  high-order  more  complex 

modes.  Preliminary  results  indicate  that  the  floor  contribution  is  neg-  | 

ligible  relative  to  that  of  the  surface  mode  for  cr„,  /a  = 10  ^ but 

floor  sea 

forms  a significant,  although  not  major,  portion  of  the  net  field  for 

a /cr  = 10  ^ . 

floor  sea 

Analysis  of  the  large  volume  of  data  generated  is  continu- 
ing. More  detailed  results  will  be  presented  in  a paper  being  prepared 
for  publication. 
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D.  Project  4504.  TROPOSPHERIC  RADIO  PROPAGATION 

Principal  Investigator:  A.  T.  Waterman 
Staff:  R.  D.  Fleming 

1.  Objective 

The  objective  of  this  research  is  to  investigate  atmospheric 
parameters  and  phenomena  in  the  lower  troposphere  based  on  radio  propa- 
gation in  the  S-  and  K-bands . The  goals  are  to 

r 

• measure  wind  velocity  in  the  common  volume  of  a transhorizon 
troposcatter  propagation  path  by  means  of  doppler  techniques 

• probe  the  structure  of  winds  and  turbulence  in  the  lower 
troposphere  under  various  conditions  of  atmospheric  stabil- 
ity and  instability 

2 . Current  Status  of  Work 

The  equipment  for  wind  measurement  on  a transhorizon  propaga- 
tion path  has  undergone  a set  of  calibration  checks  to  ensure  accuracy. 

The  plan  of  operation  is  to  involve  the  simultaneous  reception  of  3 GHz 
CW  signals,  from  a transmitter  located  100  miles  away  (beyond  line-of- 
sight),  on  each  of  two  narrow  beams  aimed  on  either  side  of  the  great- 
circle  bearing  toward  the  transmitter  by  small  amounts  (1°  or  less).  If 
a wind  component  is  transverse  to  the  path  in  the  region  of  the  common 
volume,  it  is  anticipated  that  the  signal  received  on  one  beam  will  be 
doppler  shifted  to  higher  frequencies  and  that  the  signal  on  the  other 
will  be  shifted  to  lower  frequencies.  The  magnitude  of  the  doppler  dif- 
ference should  be  proportional  to  transverse-wind  velocity.  The  two 
beams  will  be  formed  during  data  reduction  because  the  receiving  antenna 
is  a sampling  array  that  measvires  and  stores  on  digital  tape  the 
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amplitudes  and  phases  of  the  signal  as  received  on  12  individual  antennas. 
To  obtain  high  azimuthal  resolution,  the  antennas  are  arranged  in  a hori- 
zontal linear  array. 

The  linear  array  to  be  used  is  an  array  of  12  4-ft  parabolic 
dishes  mounted  on  a 70-ft  tower.  Because  the  tower  had  been  used  in  pre- 
vious experiments  in  a vertical  position,  it  was  necessary  to  move  it  to 
a horizontal  position.  The  balance  of  the  equipment  has  been  checked  and 
is  in  operating  condition. 

Recent  work  has  been  directed  toward  calibrating  the  phase  mea- 
surement system  and  organizing  the  computer-controlled  data-gathering  pro- 
grams . The  next  step  will  be  chocking  the  system  on  a llne-of-sight  path 
before  moving  to  the  transhorizon  conf igaration . 


magnetron  capable  of  producing  «1-GW  =1-Lisec  («:1  kJ  energy)  pulses  at  =3 
GHz  . 

2.  Current  Status  of  Work 

Research  conducted  elsewhere  [1,21  has  resulted  in  an  experi- 
mental cold-cathode  relativistic  magnetron  capable  of  producing  1 to  2 GW 
pulses  at  3 GHz,  with  pulse  lengths  on  the  order  of  tens  of  nanoseconds. 
During  this  reporting  period,  we  have  considered  an  improved  version  of 
this  tube,  using  a hot  cathode  and  better  vacuum  conditions,  and  incor- 
porating the  features  of  a high-power  long-anode  tube  developed  by  Boot 
[31.  Because  one  of  the  intended  applications  of  this  tube  is  to  the 
study  of  nonlinear  parametric  interactions  in  plasmas,  support  for  the 
magnetron  construction  has  been  obtained  from  the  Physics  Division  of  the 
NSF.  Our  preliminary  design  for  the  power  supply  of  the  tube  has  been 
submitted  for  technical  discussion  and  bids  to  three  manufacturers. 
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B.  Ginzton  Laboratory.  TWO-PHOTON  RESONA.NTLY  PUMPED  IR  UP-CONTERTERS 
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1.  Objective 


The  goals  of  this  project  are  the  development  and  extension  of 
efficient  III  up-conversion  techniques  in  metal  vapors,  particularly  image 
up-convers ion . 

2.  Current  Status  of  Work 

We  have  demonstrated  the  up-conversion  of  2.9  vijn  IR  images  to 
4550  A in  Cs  vapor,  A power  efficiency  of  20  percent,  with  over  1000  re- 
solvable spots,  was  achieved  using  a pump  power  of  only  8 kW  focused  in 
2 

an  area  of  0.1  cm"”.  The  Cs  cell  had  an  active  length  of  0.2  cm  and  was 

operated  at  a density  of  1.4  X 10^^/cc.  The  pumping  laser,  Nd : lanthanum- 

2 2 

berylate,  has  a natural  two-photon  coincidence  with  the  Cs  6s  S to  7s  S 
transition,  resulting  in  a simple  practical  system  with  a number  of  po- 
tential applications. 

We  have  since  designed  and  constructed  a new  Cs  IR  imaging  sys- 
tem that  has  significantly  increased  resolution  and  efficiency.  As  seen 
in  Table  6.1,  this  increase  was  accomplished  through  higher  number  den- 
sity Cs  cells  of  larger  area;  included  for  comparison  are  the  parameters 
of  a LiNbO^  up-converter  having  the  same  performance.  The  crystal  system 
requires  1250  times  more  pump  power. 

Our  most  severe  technical  problem  has  been  the  construction  of 
a reliable  long-lived  Cs  cell.  The  cell  used  in  the  demonstration  of  the 
IR  up-conversion  consisted  of  two  metalized  sapphire  windows  brazed  to 
kovar  rings.  The  rings  were  hcliarced  to  Varian  conflat  flanges  that  were 
then  bolted  together  by  means  of  standard  copper  gaskets.  A nickel  side 
arm  was  heliarced  to  the  flanges.  The  cell  appeared  to  fail  at  the  braze 
between  the  sapphire  and  kovar  ring. 

Because  we  have  encountered  difficulties  with  the  construction 
of  the  spring-loaded-window  Cs  cell  described  in  the  last  status  report, 
we  have  designed,  constructed,  and  tested  another  cell.  Incorporated  in 
this  design  are  sapphire  windows  nickel-zirconium  brazed  to  nickel  rings. 
The  disadvantage  of  these  windows  is  that  they  cannot  withstand  repeated 
thermal  cycling  because  of  the  stress  induced  by  the  differing  thermal 
expansion  coefficients  of  nickel  and  sapphire.  Outweighing  this  disad- 
vantage is  the  fact  that  the  nickel-zirconium  braze  is  extremely  resis- 
tant to  Cs  attack. 
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Table  6.1 

DESIGN  PARAMETERS  FOR  A Cs  AND  LiNbO^  2.9  ^ UP-OONVERTER 


Parameter 

Cs 

LiNbO^ 

Pump  wavelength  (p) 

1 .0790 

1.06 

IR  wavelength  (p) 

2 .94 

2.94 

Sum  wavelength  (A) 

4560 

7800 

Pump  bandwidth  (cm 

0.1 

= 1.0 

Pump  pulselength  (nsec) 

50 

50 

Photon  efficiency  (^) 

27.5 

27.5 

Resolution 

10  = 

10  = 

Pump  power  (W) 

1.44  X 10= 

9 X 10^ 

“3 

Cs  density  (cm  ) 

17 

4.44  X 10 

— 

2 

Area  (cm  ) 

3.74 

2 

Length  (cm) 

0,2 

0 .46 

The  new  windows  were  heliarced  into  Varian  conflat  flanpies  and  i 

bolted  together  using  copper  gaskets,  as  before.  A stainless-steel  base  j 

1 

was  heliarced  to  the  flanges,  and  this  entire  assembly  was  then  metal  O-  ' 

1 

ring  sealed  to  the  base  plate  of  a vacuum  chamber.  A side  arm  was  heli-  I 

arced  to  the  base  plate  and  was  also  connected  to  a hole  in  the  base  plate  | 

1 

by  means  of  a capillary  tube.  To  load  the  cell,  a clean  glass  ampoule  of 
Cs  is  inserted  in  the  side  arm.  E^vacuating  the  outer  cell  empties  the 

space  between  the  sapphire  windows  and  side  arm  via  the  capillary  tube.  j 

The  ampoule  is  then  crushed  by  flexing  the  thin-walled  side  arm.  When  ^ 

the  side  arm  is  heated,  Cs  begins  to  leak  out  through  the  capillary  tube 
but  condenses  at  the  cold  end  farthest  from  the  side  arm  and  clogs  the 

tube,  thereby  sealing  off  the  side  arm  and  cell  (Fig.  6.1).  Thermal  cy-  ^ 

cling  of  the  window  brazes  is  avoided  by  keeping  the  cell  body  continu-  ^ 

ously  hot  although  the  side  arm  is  allowed  to  cool  between  experiments.  i 
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Fig.  6.1.  1 — Cs  UP-CONVERTER  CELL. 


This  design  has  proven  reliable.  The  nickel-zirconium  brazes 

have  held  up  well,  and  there  has  been  no  evidence  of  condensation  on  the 

sapphire  windows.  The  vacuum  surrounding  the  cell  prevents  oxidation  of 

the  brazes,  copper  gaskets,  and  heater  wires. 

Using  this  cell,  we  have  increased  the  field  of  view  of  the 

up-converter  and  have  measured  resolution  of  the  images.  Figure  6.2  is 

a schematic  of  the  experimental  setup.  With  a pump  power  of  200  kW  fo- 
2 

cused  on  a 1.1  cm  area,  a f ield-of-view  of  greater  than  10,000  resolv- 
able spots  has  been  achieved.  This  is  10  times  greater  than  our  earlier 

results  and  approaches  TV  quality.  The  Cs  cell  has  an  active  length  of 

17 

0.3  cm  and  can  operate  at  a density  of  approximately  2 x 10  /cc.  The 
test  object  photographed  was  a metal  cut-out  resolution  test  target.  Tbe 
images  showed  a resolution  of  at  least  12  lines/mm,  and  the  diffraction 
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Beam  Splitter 

Teleacopc 


: 20  cm  focal  length  lena 

1*2  • 5^  focal  length  lens 

F,  : filter  for  1.06^  (i  , 1.0?9  ^ 2.9^  u 

Fig.  6.2.  EXPERIMENTAL  SETUP. 

limit  for  the  optical  system  was  approximately  28  lines/mm.  Fine  adjust- 
ment ot  the  optical  components  may  increase  image  resolution. 

Currently,  we  are  measuring  conversion  efficiency  as  a function 
of  Cs  density.  We  will  then  make  still  further  measurements  and  charac- 
terize image  quality. 
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ABSTRACT 

The  growth  of  high  resistance  liquid-phase  epitaxial  (LPE)  GaAs 
layers  requires  close-compensation  of  shallow  donors  and  acceptors  aris- 
ing from  residual  impurities  in  the  source  materials  and  from  chemical 
reactions  between  the  growth  system  components.  Once  near-compensation 
is  attained,  further  improvements  can  be  obtained  by  the  addition  of 
chromium  (Cr),  which  forms  deep  acceptor  levels,  thus  yielding  semi-in- 
sulating LPE  GaAs.  High  resistance  LPE  layers  having  sheet  resistances 

5 

of  10  ohms/'  or  more  and  low  deep  level  impurity  concentrations  can  be 
used  as  buffer  layers  in  GaAs  field  effect  transistors  (FET)  to  improve 
device  performance. 

Three  different  horizontal  tilt  LPE  growth  systems  were  investi- 
gated. These  in  order  of  investigation  were  SiO^-C-H^,  SiO^-BN(C)-H^ , 
and  SiO^-B.\-H^  systems.  It  was  found  that  the  chemical  reactions  in 
each  of  these  systems  can  be  controlled  by  systematic  bakeouts  of  the 
arsenic-saturated  gallium  melts  before  each  growth.  Low  bakeout  temper- 
atures provided  n-type  layers  while  high  bakeout  temperatures  gave  p-type 
layers  and  a critical  intermediate  bakeout  temperature,  unique  to  each 
system,  gave  high  resistivity  layers  with  close  to  1.  The  tran- 

sition  bakeout  temperatures  for  the  SiO,^-C-H^,  .S  iO,^-B.\  (C)-H,^ , and  .SiO,^- 
BN-H^  systems  were  found  to  bo  775,  700,  and  675°C,  respectively.  The 
dominant  impurities  were  found  to  be  silicon,  carbon,  and  oxygen.  Their 
relative  concentrations  were  functions  of  bakeout  temperature  and  were 
responsible  for  the  conductivity  type  transition  of  the  grown  layers. 

The  optimum  growth  conditions  for  growing  high  resistance  LPE  lay- 
ers in  each  of  these  systems  were  established.  The  optimum  compensated 
layers  grown  in  these  systems  without  deep  level  impurity  doping  had  re- 
sistivities as  high  as  700  ohm-cm  and  sheet  resistances  in  the  range  of 

5 . „ 

5 y 10  ohms . Semi-lnsulatinp  layers  with  resistivities  vip  to  o ^ 10 

9 

onm-cm  and  sheet  resistances  above  10  ohms''  were  grown  in  the  .SiO.j- 
BN(C)-H2  system  from  Ga  melts  doped  with  0.5  atm.  ' Cr  and  baked  out  at 
the  transition  temperature. 
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The  transport  related  properties  of  semi-insulating  LPE  GaAs  were 

investigated.  Photoluminescence  measurements  showed  a dominant  peak  at 

0.83  eV  in  both  Cr-doped  GaAs  commercial  substrates  and  grown  LPE  layers. 

The  temperature  variation  of  conductivity  and  the  Hall  coefficient  of  Cr- 

doped  semi-insulating  LPE  layers  were  studied  using  differential  van  der 

Pauw  techniques.  The  thermal  activation  energies  of  these  layers  were  in 

the  range  of  0.7  eV . The  compensation  ratios  in  these  layers  were  found 

to  be  in  the  range  of  1.4-2. 6,  and  the  deep  level  concentrations  were  es- 

14  -3 

timated  to  be  in  the  low  10  cm  range. 

Transient  capacitance  measurements  were  made  to  study  the  interface 

properties  of  FET  type  devices  with  and  without  LPE  buffer  layers.  No 

detectable  traps  were  found  in  the  buffered  devices  whereas  two  deep  hole 

traps,  located  0.58  and  0.81  eV  above  the  valence  band,  were  detected  in 

the  unbuffered  FET  devices  when  the  gate  depletion  region  approached  the 

substrate  interface.  The  trap  concentrations  were  1 X 10^^  and  2 X 10^^ 
-3 

cm  , respectively.  An  interface  model  was  developed  and  used  to  suc- 
cessfully explain  that  the  minority  carrier  trapping  effects  which  ap- 
peared in  the  active  layer  were  due  to  the  image  effect  of  majority  car- 
rier traps  in  the  Cr-doped  substrate  acting  through  the  interface  space- 
charge  region.  The  results  of  this  work  are  a positive  demonstration  of 
the  usefulness  of  buffer  layers  in  FET  devices  prepared  on  semi-insulat- 
ing GaAs  substrates  . 
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PERFORMANCE  BOUNDS  FOR 
PARALLEL  PROCESSORS 
by  Ruby  Bei-Loh  Lee 


JSEP  Contract 
N00014-75-C-0601 
TR  No.  125 
November  1976 


ABSTRACT 


A general  model  of  computation  on  a p-parallel  processor  is  proposed , 
distinguishing  clearly  between  the  logical  parallelism  (p*  processes) 
inherent  in  a computation  and  the  physical  parallelism  (p  processors) 
available  in  the  computer  organization.  This  shows  the  dependence  of 
performance  bounds  on  both  the  computation  being  executed  and  the  com- 
puter architecture.  We  formally  derive  necessary  and  sufficient  condi- 
tions for  the  maximum  attainable  speedup  of  a p-parallel  processor  over 


P P* 

a uniprocessor  to  be  Sp  < minCr — -)  , where  In  p approximates 


H 


th 


'In  p’ln  p* 


the  p harmonic  number.  We  also  verify  that  empirically-derived 


speedups  are  0 — -)  . Finally,  we  discuss  related  performance  measures 
of  minin 
product , 


In  p* ' 

of  minimum  execution  time,  maximum  efficiency,  and  minimum  space-time 
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ABSTRACT 

The  objective  of  this  report  is  to  describe  the  spreading  resis- 
tance method  for  determination  of  impurity  profiles,  as  implemented  in 
the  Stanford  University  Integrated  Circuit  Laboratory.  The  main  topics 
that  are  considered  include  the  measurement  technique,  analysis  of  the 

data,  accuracy  limitations  and  applications,  • 

Discussion  of  sample  preparation  focuses  on  the  relative  bevel  ]i 

quality  achieved  with  various  pwlishing  compounds.  A technique  for  ac-  , 

I 

curately  measuring  bevel  angles  on  the  order  of  16  is  outlined.  The 
non-ideai  nature  of  the  spreading  resistance  to  resistivity  calibration 
is  discussed  with  respect  to  probe  contact  characteristics,  sample  pre- 
paration and  orientation  and  impurity  type.  The  multilayer  correction 

/ 

method  and  an  efficient  implementation  algorithm  are  presented  . 

Comparison  of  spreading  resistance  with  Incremental  sheet  resis- 
tance measurements  reveals  both  the  overall  usefulness  and  accuracy  as  ' 

well  as  the  limitations  of  the  technique.  Several  causes  for  the  dis- 
crepancies in  the  comparison  are  proposed  and  future  areas  of  investiga- 
tion are  suggested.  Finally,  several  applications  are  presented,  in- 
cluding process  development  and  characterization,  device  modeling  and 
process  simulation. 


SPREADING  RESISTANCE 
FOR  IMPURITY  PROFILES 
by  Donald  C,  D'Avanzo, 
Robert  D . Rung , and 
Robert  W.  Dutton 


iii 


SKL-77-0(l  J 


I 


i 


HIGH  VOLTAGE  aiOS 
INTEGRATED  CIRCUITS 
by  Gerald  Alan  May 


Industrial  Gifts  from  Xerox 
and  NSF  Grant  ENG75-22329 
Tech  . Summary  Rep  . No  . 4834-2 
January  1977 


ABSTRACT 


I 


A new  high  voltage  silicon-gate  simultaneous  diffusion  CMOS  process 
has  been  introduced  in  this  study.  Low  voltage  (15  V)  high-density  sil- 
icon-gate CMOS  is  fabricated  on  the  same  integrated  circuit  die  as  high 
voltage  (160  V)  p-channel  transistors.  This  feature  permits  for  the 
first  time  the  design  of  high  voltage,  low  power  integrated  circuits 
needed  in  biomedical  applications  with  the  low  voltage  signal  processing 
circuitry  and  high  voltage  drive  circuitry  on  the  same  monolithic  sub- 
strate . 

Several  types  of  high  voltage  p-channel  transistors  were  developed. 
By  incorporating  a lightly  doped,  deep  p diffusion  for  the  drain  re- 
gion, a high  voltage  device  with  a thin  gate  oxide  for  high  transconduc- 
tance has  been  fabricated  to  operate  with  -160  volts  at  the  drain-sub- 
strate .junction.  In  addition,  by  utilising  a previously  reported  field 
plate  principle,  a polycrystalline  field  plate  surrounding  the  p drain 
.junction  increases  the  breakdown  of  the  high  voltage  p-channel  transis- 
tor to  -300  volts.  The  breakdown  characteristics  of  these  transistors 
have  been  explained  by  a two  dimensional  computer  aided  solution  to 
Poisson's  equation.  The  field  crowding  occurring  under  the  gate  oxide 
in  conventional  MOS  transistors  with  shallow  p^  drain  diffusions  is 
shown  to  be  reduced  in  the  deep  p diffusion.  Both  of  the  high  voltage 
transistors  can  be  easily  fabricated  using  conventional  silicon-gate 
processing  technology. 

The  high  voltage  CMOS  technology  provides  the  first  such  process  to 
make  possible  the  realization  of  all  of  the  low  power  custom  circuitry 
necessary  for  the  development  of  a compact,  light  weight,  one-hand  read- 
ing aid  for  the  blind  (Optacon).  A high  voltage  micropower  CMOS  tactile 
display  driver  multiplexer  has  been  developed  to  drive  the  piezoelectric 
transducers  in  the  tactile  screen  of  the  custom  CMOS  Optacon  system.  A 
novel  push-pull  bootstrapping  output  driver  is  employed  to  gate  100  volt 
square  waves  to  high  capacitive  loads  (such  as  piezoelectric  transducers) 
Using  Fxasitive  feedback  from  source  to  gate  of  the  high  voltage  device 
enables  the  push-pull  output  stage  to  be  driven  from  low  voltage  CMOS 
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circuitry.  Kach  output  stage  can  deliver  9 milliwatts  to  the  load  while 
consuming  less  than  500  microwatts . 

A new  trichloroethylene  (TCE)  oxidation  technology  is  incorporated 
into  the  CMOS  fabrication  process  , The  TCE  oxidation  results  in  greater 
oxide  stability  and  longer  MOS  capacitor  storage  times.  The  average  di- 
electric breakdosvn  strength  is  also  improved  and  the  frequency  of  low 
field  oxide  breakdown  is  reduced.  The  TCE  oxidation  improves  the  high 
voltage  device  properties  by  exhibiting  faster  junction  breakdown  walk- 
out with  lower  hot  electron  Injection  current  into  the  gate  oxide. 


The  first  version  of  the  Stanford  University  IC  process  simulator 
program,  SUPREM,  is  described.  The  first  part  of  this  report  outlines 
the  development  of  numerical  models  for  the  one-dimensional  redistribu- 
tion of  impurities  in  silicon  resulting  from  moving  boundaries  as  in 
oxidation  and  epitaxy.  Also  described  is  the  diffusion  model  used  in 
SUPREM  to  simulate  concentration  dependent  dif fusivities  and  cooperative 
diffusion.  The  second  part  (Appendix  1)  contains  the  operation  manual 
for  SUPREM  which  describes  in  detail  the  user-program  interfacing.  The 
last  part  of  the  report  (Appendix  2)  consists  of  an  example  of  program 
usage  to  simulate  a multi-step  process  together  with  the  resulting  out- 
puts . 


SUPPRESSION  EFFECTS 
ASSOCIATED  WITH  VLF 
TRANSMITTER  SIGNALS 
INJECTED  INTO  THE 
MAGNETOSPHERE 
by  Rajagopalan  Raghuram 


Office  of  Naval  Research 
Contract  Nonr  N00014-76-06S9 
NSF  Atmospheric  Sciences  Sec. 
Grant  ATM75-07707  Al,2 
NSF  Div . of  Polar  Programs 
Grant  DPP- 74-0493  Al,2 
TR  No.  3456-3 
March  1977 


ABSTRACT 

Magnetospheric  growth  of  coherent  VI ? signals  transmitted  from  Siple 
Station,  Antarctica  is  inhibited  by  whistler  mode  echoes  of  earlier 
transmitter  signals.  This  new  phenomenon,  called  echo-induced  suppres- 
sion, is  observed  at  least  a third  of  the  time  that  transmissions  from 
Siple  Station  are  detected  at  the  receiving  station  in  Roberval , Quebec, 
Canada.  Suppression  levels  as  high  as  20  dB  are  observed.  Though  the 
echo  is  usually  much  weaker  than  the  direct  signal,  the  level  of  suppres- 
sion is  directly  related  to  the  amplitude  of  the  echo.  The  echoes  reduce 
triggering  of  emissions  as  well  as  the  growth  of  the  signal.  Echo-induced 
suppression  is  not  explained  by  linear  wave-wave  Interference.  The  echoes 
are  thought  to  restrict  growth  by  reducing  the  coherence  of  the  total  in- 
put signal.  According  to  quasi-linear  theory,  wave  growth  is  regulated 
by  changes  in  the  particle  pitch  angle  distribution  that  results  from 
wave-particle  interaction.  These  new  results  suggest  that  coherent  waves 
tend  to  limit  their  own  growth,  even  when  modification  of  the  particle 
pitch  angle  distribution  is  unimportant.  Other  phenomena  such  as  whis- 
tler-induced suppression  of  signal  growth  are  thought  to  be  related  to 
echo-induced  suppression.  Wave-induced  growth  suppression  provides  an 
indirect  verification  of  the  existence  of  discrete  ducts  for  the  propa- 
gation of  VLF  signals  in  the  magnetosphere.  Another  suppression  phenom- 
enon, also  discovered  from  VLF  transmissions,  is  the  suppression  of  mid- 
latitude hiss  by  as  much  as  6 dB  in  a band  up  to  200  Hz  wide  iust  below 
the  transmitter  frequency.  This  quiet  band  develops  in  5 to  25  s and 
lasts  up  to  a minute  after  the  end  of  transmissions.  Tliere  i'-  evidence 
that  power-line  radiation  in  the  magnetosphere  at  harmonics  of  60  11/  also 
produces  quiet  bands.  Quiet  bands  are  explained  by  pitch  angle  scatter- 
ing of  electrons  before  they  roach  the  growth  region  near  the  equator.  A 
change  in  the  distribution  function  of  electrcins  in  the  appropriate  par- 
allel velocity  range  is  thought  to  be  responsible  for  the  (luiet  band. 
(Juiet  bands  lead  to  indirect  estimates  of  duct  widths,  scattering  regions 
and  p>erturbations  in  electron  distribution  functions. 
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SEQUENTIAL  PREFETCH 
STRATEGIES  FOR  IN- 
STRUCTIONS AND  DATA 
by  B.  Ramakrishna  Rau 


ERDA  Contract 
EY-76-S-03-0326-PA  39 
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TR  No.  131 
SU-326-P. 39-18 
January  1977 


ABSTRACT 


An  investigation  of  sequential  prefetch  as  a means  of  reducing  the 
average  access  time  is  conducted.  The  use  of  a target  instruction  buf- 
fer is  shown  to  enhance  the  performance  of  instruction  prefetch.  The 
concept  of  generalized  sequentiality  is  developed  to  enable  the  study 
of  sequentiality  in  data  streams.  Generalized  sequentiality  is  shown 
to  be  present  to  a significant  degree  in  data  streams  from  measurements 
on  representative  programs.  This  result  is  utilized  to  develop  a data 
prefetch  mechanism  which  is  found  to  be  capable  of  anticipating,  on  the 
average,  about  75'^  of  all  data  requests. 
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A NEW  MODEL  FOR  BORON 
DIFFUSION  IN  SILICON 
by  John  Ross  Anderson 


NSF  Grant  GH-31999 
TR  No.  4737-1 
May  1976 


ABSTRACT 

A new  model  for  boron  diffusion  in  silicon  is  proposed  in  which  the 
dominant  diffusing  species  is  a neutral  boron-vacancy  pair.  These  pairs, 
each  consisting  of  a negatively  charged  boron  atom  and  a positively 
charged  vacancy,  repeatedly  form  and  breakup,  establishing  a boron  dis- 
tribution in  the  silicon.  The  boron  profile  is  determined  by  the  simul- 
taneous solution  of  the  diffusion  equations  for  boron  and  boron-vacancy 
pairs,  where  the  diffusion  coefficient  of  boron  by  itself  is  assumed  to 
be  negligible.  However,  since  boron  does  get  into  the  lattice  by  the 
diffusion  and  breakup  of  boron-vacancy  pairs,  an  apparent  boron  diffu- 
sion coefficient  can  be  determined  and  is  found  to  be  equal  to  the  ratio 
of  boron-vacancy  pair  concentration  to  unpaired  boron  concentration 
times,  the  boron-vacancy  diffusivity.  It  is  this  apparent  boron  diffu- 
sion coefficient  that  appears  in  the  literature. 

The  model  also  explains  the  concentration  dependence  of  the  appar- 
ent boron  diffusion  coefficient,  although  the  boron-vacancy  diffusion 
coefficient  is  not  concentration  dependent. 

However,  the  lifetime  and  diffusion  coefficient  of  boron-vacancy 
pairs  are  not  readily  determinable  from  thermal  diffusion  data. 

The  boron-vacancy  diffusion  model  is  also  shown  to  fit  data  from 
proton-enhanced  diffusion  of  boron  in  silicon.  In  particular,  it  suc- 
cessfully explains  the  dip  in  the  electrically  active  boron  concentra- 
tion often  seen  near  the  proton  projected  range  after  a proton  bombard- 
ment. In  addition,  the  boron-vacancy  diffusion  coefficient  and  lifetime 

may  be  estimated  from  proton-enhanced  diffusion  experiments.  At  750 °C 

-14 

the  boron-vacancy  diffusion  coefficient  is  estimated  to  be  5.7  x 10 

2 , 3 

cm  /sec  and  lifetime  is  estimated  to  be  4.5  x 10  sec.  The  value  found 

for  the  boron-vacancy  diffusion  coefficient  is  about  4000  times  the  ex- 
pected value  for  the  apparent  intrinsic  diffusion  coefficient  of  boron 
at  750 °C. 

Qualitative  arguments  are  presented  to  show  that  a similar  impur- 
ity-vacancy diffusion  model  is  likely  for  donor  dopants  in  silicon.  Also 
qualitative  arguments  are  presented  showing  that  two  such  species  diffu- 
sion effects  as  emitter  dip,  emitter  push,  and  emitter  pull  can  be  de- 
scribed by  impurity-vacancy  diffusion  models. 
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A THEORETICAL  APPROACH 
TO  THE  CALCULATIONS  OF 
ANNEALED  IMPURITY  PRO- 
FILES OF  ION  IMPLANTED 
BORON  INTO  SILICON 
by  Alejandro  Chu 
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DAA-B07-75-C-13‘ll 
NSF  Grant  GH-31999 
TR  No . 4969-2 
May  1977 


ABSTRACT 

A multi-stream  diffusion  model  is  proposed  for  the  calculation  of 
the  annealing  behavior  of  boron  that  is  ion  implanted  into  silicon  at 
room  temperature  and  subsequently  annealed.  This  model  is  capable  of 
predicting  both  the  redistribution  and  the  electrical  activation  of  bor- 
on during  the  anneal,  as  a realistic  model  should.  The  calculated  re- 
sults compare  very  well  with  extensive  experimental  data  reported  in  the 
literature.  The  comparison  Includes  samples  that  are  implanted  at  room 
temperature  with  boron  in  the  dose  range  from  10^^  to  10^^  ions/cm^  and 
subsequently  annealed  in  the  temperature  range  from  800 °C  to  1000 °C. 

This  range  of  dose  and  annealing  conditions  includes  both  the  typical 
applications  of  ion  implantation  as  it  is  applied  in  the  fabrication  of 
devices  and  the  unconventional  cases  of  high  dose  implants  and  low  tem- 
perature annealing. 

The  annealing  model  is  in  essence  a diffusion  model  extended  to  in- 
clude the  following  situations: 

(1)  precipitation  of  boron,  when  the  boron  concentration 
level  e.xceeds  the  solid  solubility  limit,  in  high  dose 
implantation  cases 

(2)  trapping  of  boron  by  the  strain  fields  associated  with 
dislocation  dipoles  that  form  during  annealing  at  low 
temperature 

In  both  cases,  the  presence  of  this  immobile  and  electrically  inactive 
boron  alters  significantly  the  redistribution  and  electrical  activation 
behavior  of  boron  during  the  anneal. 

The  presence  of  boron  precipitates  makes  it  necessary  to  include 
four  species  in  the  annealing  model.  These  species  are  substitutional 
(electrically  active)  boron,  boron-vacancy  pairs  (electrically  inactive), 
positively  charged  vacancies,  and  the  aforementioned  immobile  boron.  The 
electrically  active  boron  is  assumed  to  diffuse  substitutionally  by  means 
of  random  encounters  with  neutral  vacancies.  Tlie  boron-vacancy  pair  is 
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assumed  to  diffuse  much  more  rapidly.  The  principal  interaction  among 
the  species  is  the  reaction  of  active  boron  with  positive  vacancies  to 
form  BV-pairs. 

The  parameters  in  the  model  are  a set  of  dif fusivities  and  life- 
times for  various  species  and  interaction.  These  parameters  are  esti- 
mated from  the  examination  of  each  diffusion  and  each  interaction.  The 
values  of  the  most  important  parameters  obey  simple  activation  energy 
relations.  In  other  cases,  the  modelling  was  either  not  carried  to  the 
extent  necessary  as  to  make  the  temperature  dependence  apparent,  or  else 
the  parameters  represent  the  composite  effects  of  very  complex  interac- 
tions. In  all  cases,  the  mathematical  model  and  the  parameter  set  are 
universal  in  the  sense  that  the  same  equations  and  parameters  can  be 
used  for  the  prediction  of  ordinary  diffusion,  proton  enhanced  diffu- 
sion, and  the  annealing  behavior  of  ion  implanted  boron  at  room  temper- 
ature under  conditions  that  include  precipitation  effects  and  low  an- 
nealing temperature  anomalies.  Furthermore,  in  most  cases  of  practical 
importance,  high  dose  and  low  annealing  temperature  conditions  are  of 
limited  interest  and,  as  a result,  the  three  stream  diffusion  model  with 
an  appropriate  subset  of  parameters  is  sufficient.  This  reduced  model 
is  still  very  important  and  attractive  because,  once  the  parameters  as- 
sociated with  high  dose  and/or  low  temperature  anomalies  are  deleted, 
the  remaining  parameters  have  simple  activation  energies  and  are  deter- 
mined by  the  specification  of  a single  variable;  the  temperature. 
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AT  ULTRA-LOW-FREQUENCIES 
USING  AIRBORNE  LOOP  ANTENNAS 
by  A . C.  Eraser-Smith,  D.  M. 
Bubenik,  and  O.G.Villard,  Jr 
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i . 

I 

; In  this  report  we  investigate  the  possibility  of  using  ultra-low- 


frequency  (ULF)  signals  from  airborne  loop  antennas  (i.e.,  magnetic  di- 


poles) for  air /undersea  communication.  Because  of  the  low  data  rate  at 
ULF,  communication  is  here  understood  to  mean  the  transfer  of  short  mes- 
sages of  high  Information  content. 

We  use  numerical  integration  to  calculate  the  1 Hz  total  magnetic 
field  amplitudes  in  the  sea  at  depths  from  0 m to  200  m due  to  airborne 
unit  moment  vertical  and  horizontal  magnetic  dipoles  at  altitudes  from 

100  m to  10  km.  Considering  the  magnetic  moments  m attainable  with 

, 12  9 2 

present  aircraft  power  and  payload  capability  (m  = 10  Am  - 10  Am) 

and  the  minimum  detectable  amplitude  for  a 1 Hz  sj.t,Tial  beneath  the  sea 

(~1  m/) , we  conclude  that  air/undersea  communication  at  1 Hz  is  possible 

under  the  following  illustrative  conditions:  For  a horizontal  plane  loop 

antenna  at  3000  m altitude  and  a ULF  receiver  at  a 100  m depth,  communi- 

7 2 

cation  is  possible  for  horizontal  distances  to  10  km  for  m = 10  Am 

9 2 

and  to  33  km  for  m = 10  Am.  The  coi-responding  limits  for  a vertical 
plane  loop  antenna  are  13  km  to  61  km.  It  is  also  possible,  if  desired, 
to  limit  communication  to  a comparatively  small  circular  area  directly 
beneath  the  aircraft.  Sea  floor  effects  can  alter  these  values  signifi- 
cantly, particularly  if  the  receiver  is  near  the  floor. 


A STRUCTURAL  DESIGN 
LANGUAGE  FOR  COMPUTER 
AIDED  DESIGN  OF  DIGITAL 
SYSTEMS 

by  W.  M.  vanCleemput 


TR  No.  136 
April  1977 


ABSTRACT 


In  this  report  a lang:uage  (SDL)  for  describing  structural  proper- 
ties of  digital  systems  will  be  presented.  SDL  can  be  used  at  all  levels 
of  the  design  process,  i.e.,  from  the  system  level  down  to  the  circuit 
level.  The  language  is  intended  as  a complement  to  existing  computer 
hardware  description  languages  which  emphasize  behavioral  description. 
The  language  was  motivated  partly  by  the  nature  of  the  design  process. 


TOTALLY  IMPLANTABLE  BIDIRECTIONAL 
PULSED  DOPPLER  BLOOD  FLOW  TELEM- 
ETRY: INTEGRATED  TIMER-EXCITER 

CIRCUITRY  AND  DOPPLER  FREQUENCY 
ESTIMATION 
by  Henry  V.  Allen 


PHS  Research  Grant 
POl  GM17940 
Department  of  Health, 
Education  and  Welfare 
TR  No.  4958-4 
May  1977 


ABSTRACT 

Although  blood  flow  is  one  of  the  basic  parameters  desired  in  phys- 
iological studies,  it  has  proven  to  be  elusive  to  obtain  because  of  base- 
line and  scale-factor  drift  in  currently  available  flowmeters.  A new 
totally  implantable  bidirectional  pulsed  doppler  blood  flowmeter  resolves 

these  problems.  Based  on  two  custom-designed  integrated  circuits  (ICs), 

3 

this  device  measures  2.9  x 3.8  x 0.8  cm  and  is  the  smallest  implantable 
flowmeter  realized  to  date.  This  report  describes  the  design  of  portions 
of  the  system  intended  for  long-term  implantation  in  chronic  animal  re- 
search . 

The  effects  of  noise  sources  in  the  timing  portion  of  the  flowmeter 
are  analyzed  through  a simple  first-order  approach;  this  jitter  produces 
spectral  broadening  of  doppler  returns  and  excessive  noise  near  such 
strong  reflections  as  vessel  walls.  Equations  are  developed  to  quantify 
these  effects . 

Quadrature  direction-sensing  techniques  are  compared  with  offset 
techniques.  Infinite  attenuation  of  low-frequency  artifacts  in  the  dop- 
pler returns  is  the  primary  reason  for  pursuing  the  quadrature  approach  ; 
another  major  factor  is  that  it  adds  minimally  to  the  implantable-system 
complexity  and  power  requirements . 

The  timer-exciter  IC  sets  the  ultrasonic  frequency,  repetition  rate 
of  the  ultrasonic  burst,  and  the  burst  length.  This  IC  contains  over  80 
transistors,  75  resistors,  and  300  pF  of  capacitance  for  six  integrated 
MOS  capacitors.  The  circuit  has  a power  drain  of  3.6  mA  at  2.75  V for  a 
20  kHz  repetition  rate  and  a 1.5  usee  output  burst.  Timing  jitter  is  lim- 
ited to  the  extent  that  noise  caused  by  strong  reflections  Increases  by 
only  4 dB  over  the  normal  receiver  noise,  and  spectral  broadening  is  less 
than  50  Hz.  A novel  triple  one-shot  chain  is  introduced  that  requires 
only  a single  timing  capacitor.  Rapid  resynchronization  of  the  ultra- 
sonic oscillator  provides  coherence  between  it  and  the  repetition  rate, 
and  the  capability  to  operate  either  in  the  pulsed  or  continuous-wave 
mode  allows  this  IC  to  be  used  in  both  these  types  of  doppler  ultrasonic 
flowmeters . 


A further  innovation  is  a frequency  estimator  that  employs  zero- 
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crossing  information  to  determine  the  magnitude  and  sign  of  the  doppler 
flow  signal;  unlike  other  zero-crossing  counters,  however,  this  system 
measures  the  center  frequency  of  wideband  gaussian  noise  which  is  sta- 
tistically similar  to  doppler  returns.  This  estimator  uses  four  zero 
crossings  and  is  the  first  such  processor  specifically  designed,  tested, 
and  intended  for  doppler  flow  measurements . 

Contributions  in  four  areas — timer  stability,  direction  sensing, 
integrated-circuit  design  and  development,  and  frequency  estimation — 
provide  a necessary  basis  for  producing  a miniature  pulsed  doppler  flow- 
meter. However,  these  areas  represent  only  part  of  the  total  effort  that 
has  been  required  to  accomplish  this  task.  The  works  of  R.  W.  Gill  and 
J.  W.  Knutti  of  Stanford  University  are  also  instrumental  in  the  suc- 
cessful implementation  of  this  system. 


ON  ACCURACY  IMPROVEMENT  AND 
APPLICABILITY  CONDITIONS  OF 
DIFFUSION  APPROX I^UT ION  WITH 
APPLICATIONS  TO  MODELLING  OF 
COMPUTER  SYSTEMS 
by  Philip  S . Yu 


Ballistic  Missile  Defense 
Systems  Command  Contract 
DASG60-77-C-0073 
Computer  time  at  SLAC 
TR  No.  129 
January  1977 


ABSTRACT 


Starting  with  single  server  queueing  systems,  we  find  a different 
way  to  estimate  the  diffusion  parameters.  The  boundary  condition  is  han- 
dled using  the  Feller's  elementary  return  process.  Extensive  comparisons 
by  asymptotic,  simulation  and  numerical  techniques  have  been  conducted 
to  establish  the  superiority  of  the  proposed  method  compared  with  con- 
ventional methods.  The  limitation  of  the  diffusion  approximation  is  also 
investigated.  When  the  coefficient  of  variation  of  interarrival  time  is 
larger  than  one , the  mean  queue  length  may  vary  over  a wide  range  even 
if  the  mean  and  variance  of  interarrival  time  are  kept  unchanged  . The 
diffusion  approximation  is  applicable  under  the  condition  that  the  high 
variation  of  interarrival  time  is  due  to  a large  number  of  short  inter- 
arrival times.  Case  studies  are  conducted  on  2-stage  hyperexponential 
distributions.  A similar  anomaly  is  observed  in  two  server  closed  queue- 
ing networks  when  the  service  time  of  any  server  has  a large  coefficient 
of  variation.  Again,  a similar  regularity  condition  on  the  service  time 
distribution  is  required  in  order  for  the  diffusion  approximation  to  be 
applicable.  For  general  queueing  networks,  the  problems  become  more  com- 
plicated. A simple  way  to  estimate  the  coefficient  of  variation  of  in- 
terarrival time  (when  the  network  is  decomposable)  is  proposed.  Besides 
the  anomalies  cited  before,  networks  under  certain  topologies,  such  as 
networks  with  feedback  loops,  especially  self  loops,  can  not  be  decom- 
posed into  separate  single  servers  when  the  coefficient  of  variation  of 
service  time  distributions  become  large,  even  if  the  lai’ge  variations 
are  due  to  a large  number  of  short  service  times  . Nevertheless  , the  de- 
composability  of  a network  can  be  improved  by  replacing  each  server  with 
a self  loop  by  an  equivalent  server  without  a self  loop.  Finally,  wo 
consider  the  service  center  with  a queue  dependent  service  rate  or  ar- 
rival rate.  Generalization  to  two  server  closed  queueing  networks  where 
each  server  may  have  a self  loop  is  also  considered. 
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PASSAGE  TIME  DISTRIBUTIONS 
FOR  A CLASS  OF  QUEUEING 
NETWORKS:  CLOSED,  OPEN, 

OR  MIXED,  WITH  DIFFERENT 
CLASSES  OF  CUSTOMERS  WITH 
APPLICATIONS  TO  COMPUTER 
SYSTEM  MODELING 
by  Philip  S . Yu 


Ballistic  Missile  Defense 
Systems  Command  Contract 
DASG60-77-C-0073 
Computer  time  at  SLAC 
TR  No.  135 
March  1977 


ABSTRACT 


Networks  of  queues  are  important  models  of  multiprogrammed  time- 
shared  computer  systems  and  computer  communication  networks.  Although 
equilibrium  state  probabilities  of  a broad  class  of  network  models  have 
been  derived  in  the  past,  analytic  or  approximate  solutions  for  response 
time  distributions  or  more  general  passage  time  distributions  are  still 
open  problems.  In  this  paper  we  formulate  the  passage  time  problem  as 
a "hitting  time"  or  "first  passage  time  ' problem  in  a Markov  system  and 
derive  the  analytic  solution  to  passage  time  distributions  of  closed 
queueing  networks.  Efficient  numerical  approximation  is  also  proposed. 
The  result  for  closed  queueing  networks  is  further  extended  to  obtain 
approximate  passage  time  distributions  for  open  queueing  networks.  Fi- 
nally, we  employ  the  techniques  derived  in  this  paper  to  study  the  in- 
terfault time  and  response  time  distribution  and  density  functions  of 
multiprogrammed  computer  systems.  The  effects  of  program  behavior,  de- 
gree of  multiprogramming,  size  of  main  memory,  service  time  of  paging 
devices  and  rate  of  file  I/O  requests  on  the  shape  of  distribution  func- 
tions and  density  functions  have  been  examined. 


ERDA  Contract 
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PROGRAM  BEHAVIOR  AND 
THE  PERFORMANCE  OF 
INTERLEAVED  MEMORIES 
by  B.  Ramakrishna  Rau 


One  of  the  major  factors  influencing  the  performance  of  an  inter- 
leaved memory  system  is  the  behavior  of  the  request  sequence,  but  this 
is  normally  Ignored.  This  report  examines  this  issue.  Using  trace 
driven  simulations  it  is  shown  that  the  commonly  used  assumption,  that 
all  requests  are  equally  likely  to  be  to  any  module,  is  not  valid.  The 
duality  of  memory  Interference  with  paging  is  noted  and  this  suggests 
the  use  of  the  Least-Recently-Used  Stack  Model  to  model  program  behav- 
ior. Simulation  shows  that  this  model  is  quite  successful.  An  accurate 
expression  for  the  bandwidth  is  derived  based  upon  this  model. 
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PERFORMANCE  ANALYSIS  OF 
COMPUTER  COMMUNICATION 
NETWORKS  VIA  RANDOM 
ACCESS  CHANNELS 
by  Kiilip  S.  Yu 


Ballistic  Missile  Defense 
Systems  Command  Contract 
DASG60-77-C-0073 
Computer  time  at  SLAC 
TR  No.  137 
April  1977 


ABSTRACT 


The  field  of  computer  communication  networks  has  grown  very  rapidly 

in  the  past  few  years.  One  way  to  communicate  is  via  multiple  access 

broadcast  channels  . A new  class  of  random  access  schemes  referred  to  as 

the  M -persistent  CSMA  scheme  is  proposed.  It  incorporates  the  nonper- 
P 

sistent  CSMA  scheme  and  l-p)ersistent  CSMA  scheme,  both  slotted  and  un- 
slotted versions,  as  its  special  cases  with  p = 0 and  1,  respectively. 

The  performance  of  the  M^-pers istent  CSMA  scheme  under  packet  switching 
is  analyzed  and  compared  with  other  random  access  schemes.  By  dynami- 
cally adjusting  p,  the  unslotted  version  can  achieve  better  performance 
in  both  throughput  and  delay  than  the  currently  available  unslotted  CSMA 
schemes  under  packet  switching.  Furthermore,  the  performance  of  various 
random  access  schemes  under  message  switching  is  analyzed  and  compared 
with  that  under  packet  switching.  In  both  slotted  and  unslotted  versions 
of  the  Mj^-pers istent  CSMA  scheme,  the  performance  under  message  switching 
is  superior  to  that  under  packet  switching  in  the  sense  that  not  only  the 
channel  capacity  is  larger  but  also  the  average  number  of  retransmissions 
per  successful  message  under  message  switching  is  smaller  than  that  per 
successful  packet  under  packet  switching.  In  dynamic  reservation  schemes, 
message  switching  leads  to  larger  channel  capacity.  However,  in  both 
slotted  and  unslotted  versions  of  the  ALOHA  scheme,  the  channel  capacity 
is  reduced  when  message  switching  is  used  instead  of  packet  switching. 

This  phenomenon  may  also  happen  in  the  M -persistent  CSMA  scheme  as  p 

P 

deviates  from  0 to  1 for  certain  distributions  of  message  length.  Hence, 
the  performance  under  message  switching  may  be  superior  to  or  inferior  to 
that  under  packet  switching  depending  upon  the  random  access  scheme  being 
used  and  the  distribution  of  message  length  (usually  a large  coefficient 
of  variation  of  message  length  implies  a large  degradation  of  channel  ca- 
pacity in  this  case)  for  certain  random  access  schemes.  Nevertheless,  for 
radio  channels,  message  switching  can  achieve  larger  channel  capacity  if 
appropriate  CSMA  schemes  are  used  . A mixed  strategy  which  is  a combina- 
tion of  message  switching  and  packet  switching  is  proposed  to  improve  tl>e 
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performance  of  a point  to  point  computer  communication  network  when  its 
terminal  access  networks  communicate  via  highly  utilized  radio  channels. 
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PROPERTIES  AND  APPLICATIONS 
OF  THE  LEAST-RECENTLY-USED 
STACK  MODEL 
by  B.  Ramakrishna  Rau 


ERDA  Contract 
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ABSTRACT 


The  Least-Recently-Used  Stack  Model  (LRUSM)  is  known  to  be  a good 
model  of  temporal  locality.  Yet,  little  analysis  of  this  model  has  been 
performed  and  documented.  Certain  properties  of  the  LRUSM  are  developed 
here.  In  particular,  the  concept  of  the  Stack  Working  Set  is  introduced 
and  expressions  are  derived  for  the  forward  recurrence  time  to  the  next 
reference  to  a page,  for  the  time  that  a page  spends  in  a cache  of  a 
given  size  and  for  the  time  from  last  reference  to  the  page  being  re- 
placed. The  fault  stream  out  of  a cache  memory  is  modelled  and  it  is 
shown  how  this  can  be  used  to  partially  analyze  a multilevel  memory 
hierarchy.  In  addition,  the  Set  Associative  Buffer  is  analyzed  and  a 
necessary  and  sufficient  condition  for  the  optimality  of  the  LRU  re- 
placement algorithm  is  advanced  . 
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ATMOSPHERIC  PROBING 
USING  WAVE  PROPAGA- 
TION MEASUREMENTS 
by  A.  T.  Waterman,  Jr. 
and  R.  W.  Lee 


ABSTRACT 

The  objective  of  this  research  was  the  development  of  a technique 
for  the  remote  measurement  of  wind  and  atmospheric  turbulence.  The  tech- 
nique involves  the  transmi.ss ion  of  millimeter-length  radio  waves  or  of 
optical  waves  (or  even  of  a"oustlc  waves)  between  two  locations  in  the 
atmosphere.  By  special  modification  and  processing  of  the  transmitted 
waves,  the  measurement  is  made  sensitive  to  atmospheric  conditions  at 
iust  one  point  on  the  transmission  path,  and  to  I’ust  one  scale  size  of 
atmospheric  turbulence.  Altering  this  modification  of  the  waves  leads 
to  an  alteration  in  the  location  of  the  sensitive  point  along  the  path, 
or  to  an  alteration  in  the  sensitive  scale  size.  In  this  way,  an  exam- 
ination of  conditions  along  the  path  can  be  made;  for  example,  one  could 
measure  a profile  of  the  wind  (actually  the  transverse  component  of  the 
wind)  along  a line  between  two  spaced  points.  The  technique  falls  in  the 
category  of  remote  probing  since  no  instruments  are  placed  at  the  point 
where  the  measurement  is  made.  As  an  example,  the  wind  over  an  airport 
runway  can  be  measured  by  placing  the  transmitting  and  receiving  points 
well  to  either  side  of  the  runway. 

The  cru.x  of  the  technique  lies  in  the  particular  modification  and 
processing  of  the  transmitted  waves.  It  is  the  use  of  spatial  filters. 

In  their  crudest  form,  spatial  filters  may  be  no  more  than  gratings.  Two 
variants  of  the  technique  were  made  the  subject  of  this  research.  In  one, 
a spatial  filter  was  used  at  the  receiving  end  of  the  path  only,  and  was 
supplemented  by  a variable  temporal  filter.  In  the  second,  spatial  fil- 
ters were  used  at  both  transmitting  and  receiving  ends  of  the  path,  llio 
first  of  these  was  applied  to  a path  28  km  in  length  using  25  GHz  radio 
waves  (of  8.6  millimeter  wavelength),  since  equipment  which  could  be  used 
for  this  purpose  had  been  constructed  and  set  up  under  a previous  research 
program.  The  second  was  applied  to  shorter  paths,  ranging  from  10  to  225 
meters,  and  used  Incoherent  optical  waves  (visible  light  from  an  incan- 
descent bulb ) . 

Successful  operation  of  the  technique  was  clearly  demonstrated.  Wind 
measurements  ojtained  using  the  optical,  double-spatial-filter  approach 
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were  chocked  against  anemometer  readings.  The  longer-path-length , mil- 
limeter-wave approach  also  was  successful,  though  less  easy  to  verify. 
Attempts  to  operate  the  optical  equipment  on  a vertical  path  on  San  Fran- 
cisco's Mt . Sutro  Tower  were  less  than  satisfactory,  being  hampered  by 
intense  interference  from  the  commercial  TV  and  FM  transmitters  on  the 
tower.  Other  proposed  extensions  of  the  technique  could  not  be  pursued 
owing  to  funding  limitations.  However,  the  experimental  and  measurements 
program  was  supplemented  by  a rather  thorough  theoretical  analysis  of  the 
technique  . 
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INTERPRETIVE  MACHINES  JSEP  Contract 

by  John  K.  Iliffe  N00014-75-C-0601 

Course  at  Swansea 
TR  No.  149 
June  1977 


ABSTRACT 

These  lectures  survey  attempts  to  apply  computers  directly  to  high 
level  languages  using  microprogrammed  interpreters.  The  motivation  for 
such  work  is  to  achieve  language  implementations  that  are  more  effective 
in  some  measure  of  translation,  execution  or  response  to  the  user  than 
would  otherwise  be  obtained.  The  implied  comparison  is  with  the  estab- 
lished technique  of  compiling  into  a fixed  general-purpose  machine  code 
prior  to  execution.  It  is  argued  that  while  substantial  benefits  can  be 
e^'pected  from  microprogramming  it  does  not  represent  the  best  approach 
to  design  when  the  contributing  factors  are  analysed  in  a general  system 
context,  that  is  to  say  when  wide  performance  range,  multiple  source 
language,  and  stringent  security  requirements  have  to  be  satisfied.  An 
alternative  is  suggested,  using  a combination  of  interpretation  and' a 
primitive  instruction  set  and  providing  security  at  the  microprogram 
level . 


COMPUTING  IN  STORE 
by  John  K.  Iliffe 


JSEP  C ".tract 
N00014-75-C-0601 
TN  No,  117 
June  1977 


ABSTRACT 


These  notes  provide  an  introduction  to  the  class  of  single-instruc- 
tion, multiple-data  stream  computers  with  the  simplest  processing  ele- 
ments. Design  principles  are  explained  in  terms  of  hypothetical  Dis- 
tributed Processor  Arrays,  with  examples  drawn  from  experimental  systems. 
Emphasis  is  placed  on  (a)  minimizing  the  cost  differential  when  the  DPA 
is  compared  with  conventional  main  storage,  and  (b)  designing  the  array 
control  unit  to  support  advanced  forms  of  protection  and  language  imple- 
mentation. The  influence  of  the  DPA  on  general  system  design  is  examined 
briefly. 
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BROADCAST  PROTOOOI5 
IN  PACKET  SWITCHED 
COMPUTER  NETWORKS 
by  Yogen  K.  Dalai 
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ABSTRACT 

This  thesis  investigates  the  design  and  analysis  of  broadcast  rout- 
ing algorithms  for  use  in  store-and-forward  packet  switched  computer 
networks.  Broadcast  routing  is  taken  here  to  be  a special  case  of  multi- 
destination routing,  in  which  a packet  is  delivered  to  all  destinations 
rather  than  to  some  subset . 

We  examine  five  alternatives  to  transmitting  separately  addressed 
packets  from  the  source  to  the  destinations.  The  algorithms  are  compared 
qualitatively  in  terms  of  memory  requirements,  ease  of  implementation, 
adaptiveness  to  changing  network  conditions,  and  reliability.  The  algo- 
rithms are  also  compared  quantitatively  in  terms  of  the  number  of  packet 
copies  generated  to  perform  broadcast  and  the  delays  to  propagate  the 
packet  to  all  destinations . Lower  bounds  on  the  performance  measures 
are  determined  for  all  the  algorithms  by  examining  regular  graphs. 

Protocols  that  provide  reliable  communication  using  broadcast  rout- 
ing, 1 .e . , broadcast  protocols,  are  analogous  to  Interprocess  communica- 
tion protocols  except  that  communication  is  between  one  process  and  many 
processes.  Reliable  broadcast  protocol  design  is  faced  with  problems 
similar  to  those  in  the  design  of  interprocess  communication  protocols  - 
addressing,  sequencing,  duplicate  detection  and  guarantee  of  delivery. 
This  area  presents  many  subjects  for  future  research. 

We  describe  a few  applications  for  broadcast  protocols  in  distrib- 
uted computing  environments.  In  particular,  we  show  in  detail  how  the 
catalog  of  a distributed  file  system  could  be  structured  in  a simple  way, 
if  the  system  could  make  use  of  efficient  reliable  broadcast  protocols. 

The  properties  of  reliable  broadcast  protocols  at  the  host  level 
emerge  from  the  reliability  of  the  routing  algorithms  and  the  applica- 
tions for  the  protocols.  We  have  examined  the  tradeoffs  between  global 
and  subgroup  broadcast  routing.  One  conclusion  we  offer  is  that  commu- 
nication subnets  should  support  both  capabilities  in  the  form  of  multi- 
destination  addressing  and  reverse  path  forwarding,  respectively. 

An  outcome  of  the  investigation  of  broadcast  routing  algorithms  is 
the  formulation  of  two  distributed  (parallel)  algorithms  for  constructing 


minimal  spanning  trees.  We  believe  that  these  algorithms  are  the  first 
of  their  kind.  The  formulation  of  such  algorithms  has  made  the  problems 
affecting  the  design  of  distributed  algorithms  in  network  environments 
clearer.  These  minimal  spanning  tree  algorithms  can  be  used  in  broad- 
cast routing,  as  well  as  other  networks  like  the  Packet  Radio  Network. 
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SYSTEM  FOR  PRINTED  CIRCUIT 
BOARD  DESIGN  USER'S  GUIDE 
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ERDA  Contract 
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June  1977 


ABSTRACT 


The  SPRINT  system  for  the  design  of  printed  circuit  boards  Is  a col- 
lection of  programs  that  allows  designers  to  Interactively  design  two- 
sided  boards  using  a Tektronix  4013  graphics  terminal.  The  major  parts 
of  the  system  are:  a compiler  for  SDL,  the  Structural  Design  Language, 
an  Interactive  component  placement  program,  an  Interactive  manual  conduc- 
tor routing  program,  an  automatic  batch  router,  a via  elimination  program 
and  a set  of  artwork  generation  programs . 
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OXIDATION  AND  EPITAXY 
by  R.  W.  Dutton,  D.  A. 
Antonladls,  J.  D.  Melndl , 
T . I . Kamlns , K . C . 
Saraswat,  B.  E.  Deal, 
and  J.  D.  Plummer 


ABSTRACT 


The  first  order  process  models  for  silicon  epitaxy  and  oxidation 
are  described.  Epitaxial  dopant  inclusion,  autodoping  and  transient 
effects  are  discussed,  and  experimental  results  are  presented.  Silicon 
orientation,  surface  doping,  and  ambient  effects  are  considered  for 
silicon-oxidation  rates . 
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ABSTRACT 


The  Gunmel-Poon  model  is  reviewed  with  emphasis  directed  toward 
automated  measurement.  An  algorithm  for  model  paramet»^r  extraction  is 
given  and  measured  and  simulated  results  are  compared . A statistical 
sample  of  area-scaled  transistors  across  a wafer  has  been  characterized. 
Factor  Analysis  is  used  to  help  reduce  the  data  and  define  a statistical 
model.  Computer  programs  for  Process  and  Device  Simulation  are  used  to 
facilitate  understanding  of  the  observed  statistical  variations  and  pa- 
rameter correlations. 
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